{ "cells": [ { "cell_type": "markdown", "id": "ca6faffd-32ed-4c32-9921-e424c7a01464", "metadata": {}, "source": [ "Magnetic field interpolation (J-TEXT)\n", "=====================================" ] }, { "cell_type": "code", "execution_count": 1, "id": "7aea8032-2ea7-4cdd-9321-ddca1e3592de", "metadata": {}, "outputs": [], "source": [ "import fusionsc as fsc\n", "from fusionsc.devices import jtext\n", "import numpy as np\n", "from tqdm.auto import tqdm\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "id": "f5e454a6-c699-475a-aa7d-d57f0fa034c8", "metadata": {}, "source": [ "Here we show a simple example to evaluate the magnetic field at desired points. We use the function \"MagneticConfig.interpolateXyz\" to evaluate our field." ] }, { "cell_type": "code", "execution_count": 2, "id": "b1219848-1e53-4923-b13b-ef417d749637", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Help on AsyncMethodDescriptor in module fusionsc.magnetics:\n", "\n", "interpolateXyz(self, points, grid=None) -> Any\n", " Evaluates the magnetic field at target positions by first computing the magnetic field\n", " at the target points (if not yet done), and then \n", " \n", " Parameters:\n", " - points: A numpy-array of shape [3, ...] (at least 1D) with the points in x, y, z coordinates.\n", " - grid: An optional grid parameter required if the field is not yet computed. The grid\n", " \n", " Returns:\n", " A numpy array of shape points.shape with the field as x, y, z field (cartesian).\n", " \n", " *Note* Has :ref:`asynchronous variant` '.asnc(...)' that returns Promise[...]\n", "\n" ] } ], "source": [ "help(fsc.magnetics.MagneticConfig.interpolateXyz)" ] }, { "cell_type": "code", "execution_count": 3, "id": "222191e1-bb47-4c83-8bb3-bff1f24ffdb6", "metadata": {}, "outputs": [], "source": [ "efitExample = jtext.exampleGeqdsk()\n", "field = fsc.magnetics.MagneticConfig.fromEFit(efitExample).compute(jtext.defaultGrid())" ] }, { "cell_type": "code", "execution_count": 4, "id": "bfa4a0a1-82cf-46cc-9ad1-023ffe11febb", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjUAAAGdCAYAAADqsoKGAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAA9hAAAPYQGoP6dpAACBwklEQVR4nO3dd1hT59sH8G8SRlhhTwUVF05EVMQ9UBx111336rCtYmulb+uotmr1Z61Va22r1rbWVbVOHFjFgQMUJw4UZcheYWae9w80kBLmySK5P73OJeec55zcSSG580wOwzAMCCGEEELqOa6uAyCEEEIIUQdKagghhBBiECipIYQQQohBoKSGEEIIIQaBkhpCCCGEGARKagghhBBiECipIYQQQohBoKSGEEIIIQbBRNcB6IJcLserV69gY2MDDoej63AIIYToMYZhkJ+fDw8PD3C5mqsLKCkpgVgsZn0fMzMz8Pl8NURU/xhlUvPq1St4enrqOgxCCCH1SGJiIho2bKiRe5eUlKBJI2ukpstY38vNzQ3x8fFGmdgYZVJjY2MDoPQXVCAQ6DgaQggh+kwoFMLT01Px2aEJYrEYqekyxEc3gsCm7rVBwnw5mvi/hFgsrlFSExERgXXr1iE6OhopKSk4fPgwRo4cWeU1f/75J7799ls8ffoUtra2GDx4MNatWwdHR8c6x60uRpnUvGlyEggElNQQQgipEW10VxDYcFklNbVVWFgIX19fzJw5E6NHj662/JUrVzB16lR89913GDZsGJKTk/Huu+9izpw5OHTokBYirppRJjWEEEKIPpIxcshYLDMtY+S1Kj948GAMHjy4xuUjIyPRuHFjfPTRRwCAJk2aYN68eVi7dm2tHldTaPQTIYQQoifkYFhvQGmTWflNJBKpJb7AwEAkJibi5MmTYBgGaWlpOHjwIIYMGaKW+7NFSQ0hhBCiJ+Rq+A8APD09YWtrq9hWr16tlvi6d++OP//8E+PHj4eZmRnc3Nxga2uLLVu2qOX+bFFSQwghhBiYxMRE5OXlKbbQ0FC13Pfhw4f4+OOPsXTpUkRHRyMsLAwvXrzAu+++q5b7s0V9agghhBA9IWMYyJi6d6p5c62mBsKsXr0a3bt3x6effgoAaN++PaysrNCzZ0+sWrUK7u7uan/M2qCkhhBCCNET5fvF1PV6TSoqKoKJiXLqwOPxAJROUqhrlNQQogNSuRximRQimQwi6et/ZVKIZFI0sBbAwcJS1yESQoxAQUEB4uLiFPvx8fGIiYmBg4MDvLy8EBoaiuTkZOzevRsAMGzYMMyZMwc//vgjgoODkZKSggULFqBLly7w8PDQ1dNQoKSGEDVjGAbfR0ci/OWz/yQsbxIYqcrqZR6HgwWduuE9vwAdRE0I0QdyMJBpsaYmKioKffv2VeyHhIQAAKZNm4Zdu3YhJSUFCQkJivPTp09Hfn4+Nm/ejEWLFsHOzg79+vXTmyHdHEYf6ou0TCgUwtbWFnl5eTT5HtGIYokEc8KO4HLyyxqVb2Jrj+/6DUEHV922RxNCKtLGZ8abx3j2yA02LCbfy8+Xo6lPqtF+vtHoJ0I0wMLUFMt79IM936LashNbtceJt6dQQkMIISxR8xMhapScL8TJ549x4tkTxKSnVFnWgW+BtX2CMaBxMy1FRwjRd+oa/WSsKKkhhKXkfCFOPX+C488eV5vIvNHHswm+7TsILpZWGo6OEFKfyF9vbK43ZlppftqyZQsaN24MPp+PgIAA3Lhxo9KyP//8M3r27Al7e3vY29sjKCioQnmGYbB06VK4u7vDwsICQUFBePr0qaafBiEVMAyD8Uf3YlXkhQoJjRmXhwGNm6GBdVm7tjnPBF/16I+dQ0ZTQkMIIWqm8aRm3759CAkJwbJly3Dr1i34+voiODgY6enpKstfuHABEydOxL///ovIyEh4enpi4MCBSE5OVpT59ttvsWnTJmzbtg3Xr1+HlZUVgoODUVJSoumnQ4gSDoeDId4tFPtmXB6CGjfFxn5DED39ffw8aCS4r1f2be3oguNj3sHUtn5aWe2XEFL/yF6PfmKzGTONj34KCAhA586dsXnzZgCAXC6Hp6cnPvzwQyxZsqTa62UyGezt7bF582ZMnToVDMPAw8MDixYtwieffAIAyMvLg6urK3bt2oUJEyZUe08a/UTU6W56KjZFR2Jo05bo36gpBObminNimQytf/0es9r7Y1HnHjB7PUkVIaT+0Obop7sPXViPfmrfOt1oP9802qdGLBYjOjpaac0JLpeLoKAgREZG1ugeRUVFkEgkcHBwAFA6MVBqaiqCgoIUZWxtbREQEIDIyEiVSY1IJFJaoVQoFNb1KRFSQXsXN/wyeJTKc/liEX4f+jYCG3hpOSpCSH1EfWrY0WjzU2ZmJmQyGVxdXZWOu7q6IjU1tUb3+Oyzz+Dh4aFIYt5cV5t7rl69Wmm1Uk9Pz9o+FULqxNHCkhIaohYlxWKIxVKlY7H3EnUUDSH6Sa/nqVmzZg327t2Lw4cPg8/n1/k+oaGhSquVJibSGwEhpH7ZtPo4Fkz/BQkvMhTHtq47hasXYnUYFVE3OTiQsdjkMO7+ehpNapycnMDj8ZCWlqZ0PC0tDW5ublVeu379eqxZswZnzpxB+/btFcffXFebe5qbmytWLNXUyqWEEKIp50/dRfjJu3j2JBXz39mOU0eiIZfLkRCfgXXLjuBVUrauQyRqImfYb8ZMo0mNmZkZ/P39ER4erjgml8sRHh6OwMDASq/79ttvsXLlSoSFhaFTp05K55o0aQI3NzelewqFQly/fr3KexJCSH0lk8lhbl7aBVJUIsHGVcfwxUd/oqRYgqJCEVYt3g9RiUTHURKiexpvfgoJCcHPP/+M3377DbGxsXjvvfdQWFiIGTNmAACmTp2q1JF47dq1+PLLL7Fjxw40btwYqampSE1NRUFBAYDSIbQLFizAqlWrcPToUdy7dw9Tp06Fh4cHRo4cqemnQwghWjfgrQ744fe5aNK8rC9h9LVnip+fPUnF1vWndBEaUTM2TU9vNmOm8RmFx48fj4yMDCxduhSpqano0KEDwsLCFB19ExISwOWW5VY//vgjxGIx3n77baX7LFu2DMuXLwcALF68GIWFhZg7dy5yc3PRo0cPhIWFsep3Qwgh+qyRtws27ZqNX384hyN7r1c4H3bkFtr4emLgMD8dREfUhW1iYuxJDa3STf1rCCH1RMKLDOze9i8unXuo8ry5uQk27poN7+ZV91kktaPNeWquPnCHNYt5agry5ejWJsVoP9/0evQTIYSQUtciHuPDKdsrTWgAQCSSYtVn+1FYQLOr11dyhsN6M2aU1BBCSD3QtVdL7Dm1CJ+uGIXO3ZqBy1P94ZWckI3vVh6FEVbCGwTqU8MOrdJNCCH1hJU1H0FDfRE01Bd5uYW4dO4hLpy5j/u3X6J8DnMp/CGO7L2OURO76i5YQnSAkhpCCKmHbO2s8NbbnfHW252RkZaHi2cf4MLpe3gaW7pa/M8bz6BlmwZo3Z5mUK9PZOBCxqIRRabGWOojSmoIIaSec3a1xdvvdMPb73RDckIWLpy5jwun7+Hr0APY8sc82Nlb6TpEUkMMy34xjJH3qaGkhhBCDEgDL0dMnt0bk2b1wvOnaXj+JBUdA5rqOixSQzSkmx1KagghxABxOBw0bUFDu4lxoaSGEEII0RMyhgsZw6JPjZEPeqOkhhBCCNETcnAgZ9FRWA7jzmponhpCCCGEGASqqSGEEEL0BHUUZoeSGkIIIURPsO9TQ81PhBBCCCH1HtXUEEIIIXqitKNw3ZuQ2FxrCCipIYQQQvSEnOUyCTT6iRBCCCHEAFBNDSGEEKInqKMwO5TUEEIIIXpCDi5NvscCJTWEEEKInpAxHMhYrLTN5lpDQH1qCCGEEGIQqKaGEEII0RMylqOfZNT8RAghhBB9IGe4kLPoKCw38o7C1PxECCGEEINANTWEEEKInqDmJ3YoqSGEEEL0hBzsRjDJ1RdKvUTNT4QQQggxCFRTQwghhOgJ9pPvGXddhXE/e0JIvRTzKAm/Ho7Es8RMMNWM9khKy4Vcbtz9DEj98WaZBDabMaOaGkJIvXM84gGORzzAz39HoqGrHXp3aoY+nZqhTVN3cLnK/RGu3X2Bi9Fx+GJOMFwdbXQUMSFEG4w7pSOE1DtSmRyXbj1T7Cel5eLPE1GYs2Ivhn+0Hd/uDMf1ey8hlcoAAI09HHDzfgImLfkNJyIeVFuzQ4guycFhvRkzSmoIIfVO6KwBGNyjFWwszZWOZ+YW4lD4HXy89m8Mfn8blv94Ci9eZQMACovFWLn9ND7beBRZeYW6CJuQamm7+SkiIgLDhg2Dh4cHOBwOjhw5UmX56dOng8PhVNjatGnD4lmrDzU/EULqFRMeF306N0efzs0hlcpw61ESLtyMQ0R0HDJzy5KV/CIRwq7EIuxKrNL1EdHPcPfJKyye0R/9urTQdviEVIn9PDW1u7awsBC+vr6YOXMmRo8eXW3577//HmvWrFHsS6VS+Pr6YuzYsbWOVRMoqSGE1FsmJjx0adsIXdo2wifT+uHh8xT8ezMOF6PikJSWW+l1ufnF+HzTcQR388Giaf0gsOJrL2hC9MjgwYMxePDgGpe3tbWFra2tYv/IkSPIycnBjBkzNBFerWml+WnLli1o3Lgx+Hw+AgICcOPGjUrLPnjwAGPGjEHjxo3B4XCwcePGCmWWL19eoerLx8dHg8+AEKLvuFwO2jbzwIcTe+HA+hn4c/VUTB3WucprTl99hMlLfsO1uy+0EyQh1ZAzHNYbAAiFQqVNJBJpJN5ff/0VQUFBaNSokUbuX1saT2r27duHkJAQLFu2DLdu3YKvry+Cg4ORnp6usnxRURG8vb2xZs0auLm5VXrfNm3aICUlRbFdvnxZU0+BEFLPcDgcNHC1xcPnadWWzcgpxIJvD2HtjnMoKhFrITpCKid/3fxU1+3NPDWenp6KWhVbW1usXr1a7bG+evUKp06dwuzZs9V+77rSePPThg0bMGfOHEXV1LZt23DixAns2LEDS5YsqVC+c+fO6Ny59NuVqvNvmJiYVJn0EEKMl1giRej3xxD1IKHG1xw+fxc37r/El3OD0cGnoQajI0TzEhMTIRAIFPvm5uZVlK6b3377DXZ2dhg5cqTa711XGk1qxGIxoqOjERoaqjjG5XIRFBSEyMhIVvd++vQpPDw8wOfzERgYiNWrV8PLy0tlWZFIpFT1JhQKWT02IUS/cblcfPX+EMjkDGQyOWQyOaTy0n8VxxT7ckhlZcdEEikkUhlMTXi6fhrECMkZLuQsJtB7c61AIFBKatSNYRjs2LEDU6ZMgZmZmcYep7Y0mtRkZmZCJpPB1dVV6birqysePXpU5/sGBARg165daNmyJVJSUrBixQr07NkT9+/fh41Nxcm1Vq9ejRUrVtT58Qgh9YsJjwsb6vxL6iEZOJCxmGuGzbW1cfHiRcTFxWHWrFlaebyaqpfz1AwePBhjx45F+/btERwcjJMnTyI3Nxf79+9XWT40NBR5eXmKLTExUcsRE0IIIfqnoKAAMTExiImJAQDEx8cjJiYGCQmlTbehoaGYOnVqhet+/fVXBAQEoG3bttoMt1oaralxcnICj8dDWppyZ720tDS19oexs7NDixYtEBcXp/K8ubm5RtoTCSGEEHVSV/NTTUVFRaFv376K/ZCQEADAtGnTsGvXLqSkpCgSnDfy8vLw999/4/vvv69znJqi0aTGzMwM/v7+CA8PV3QkksvlCA8Px/z589X2OAUFBXj27BmmTJmitnsSQggh2iYDuyYkWS3L9+nTp8qlQ3bt2lXhmK2tLYqKimr5SNqh8dFPISEhmDZtGjp16oQuXbpg48aNKCwsVIyGmjp1Kho0aKAYbiYWi/Hw4UPFz8nJyYiJiYG1tTWaNWsGAPjkk08wbNgwNGrUCK9evcKyZcvA4/EwceJETT8dQgghhOgpjSc148ePR0ZGBpYuXYrU1FR06NABYWFhis7DCQkJ4HLLqstevXoFPz8/xf769euxfv169O7dGxcuXAAAJCUlYeLEicjKyoKzszN69OiBa9euwdnZWdNPhxgBGSMHF6WTOhJCiDZpu/nJ0HAYI1yyVigUwtbWFnl5eRod8kbqF6lcjuNJ95BQkI2PWvet/gJCiFHQxmfGm8cIjRwEvrVpne9TUiDB6sAwo/18o7WfiNGTyuU4kXQPPz6KQGqxECcH1K6/F8MwVKtDCFELBhzIWfSpYbQ0pFtfUVJDjJaMkeNk0n1sfRSBFwVZAID3fXrBw9K2mitLPcxJw/7nMZjfpgec+FaaDJUQQkgNUFJDjI6MkeNU0gNsfXQR8a+TGQBw5dtgVvPuVV4rZxj8+yoOOx7fwLX0l1jhH0wJDSFEbWQMFzIW/WLYXGsIKKkhRkPGyBGW9ABbH0XgeUFmhfOftB0ASxPV030XScU4FH8PO5/cxIv8bABAR6eGmNSso0ZjJoQYl/Irbdf1emNGSQ0xeHKGQVhyac3Ms/yKyQwAdHBoiKENK86MmVqUj9+fRuGvZ7eRJy5RHDfj8rC68xBwqS8NIbVCfdCIJlFSQwzaxdSnWHf/LJ7lZ1RZ7vP2g5TeaO9lp2Dn4xs4kRALKSOvUP791t3RzNZJ7fESYshePn6F7z7ejUmLhqLLgHaK43K5XGlqD2MmAxcyFisYsbnWEFBSQwxaF6fGeK9lT/yTeBdX0p5BjoozGIz08kU7+waQyeUIf/UUOx/fwI2MytcHa27rhHmtAjUZNiEG5/rpu1g1YxskYik2LfoDP11ZDisbC0SF3wfXhIuOvVvrOkS9QM1P7Bh3SkcMnoWJKYZ6tsO4xv7gqfgmaMkzxcLW/QEAV9Je4Lt7EVUmNBwAqzsPhRmPp6mQNUJYUoJ7aWmIz8nRdSjESLUOaAprO0sAQOarHOz86hAkYim2/d8+3DxzT8fREUNBNTXE4B1NuIvPbx2BTMU8k3Nb9oSLhQ0AoJe7N3q6NcHl1Hi8e/kgSmTSCuWntugEP6cGGo+5tmRyOVILCvAyNxeJuXlIyMtFQm4eEvNKf84rEWFQ8+b43+BBug6VGCkbOyvMXzcZK6f9CAA4vvMi8nMKkRSXBkbOYN7X43UcoX6Qgws5i/oGNtcaAkpqiEHbGx+Fr2JOKBqdmto4oYm1E86lPEIDSztMb6bcjJQjLsYPDy6rTGg8LAVY1K6P5oOuxsP0dFxNSEBCueTllVAIibxi3583Zvn7Y0nvXtSxmehU96F+6DnCH5f+iQYAXDwSBQBIfp6O5GdpaNDUVZfh6QUZw4GMRRMSm2sNASU1xGD9+uQK1j84p9hvZeuGX7q/g7s5yTiX8giL2w6AOa/sTyBemIWZEfuRUFDWRGPC4So6Cq/sNBhWpqqHfGuTp60tbiYl49yzZ9WW5XI4+LJvH0wtt54aIdomzC5AVPgDZLzKrrTMjbP3MIqSGsKScddTEYPEMAy+f3heKaHp6OCJ33pOg4O5Fbo6N0Ev1+YY4NFKcf56+kuMOfebIqHhcjhY7j8Qb3mVdl4c0agN+ng01e4TqYSNuTl+HDEc0ztWnajwTUzw4/DhlNAQnbO2s0RiXCp2rjysqKX5r5vnqF8NUNZRmM1mzKimhhgUOcNgzb0w/P7shuJYd5em2BQwTjGxHp9niv91HqMYwn04/h5Cb55QNN9YmZjh+24j0dejGcy5MYhIfYYvOg7Q/pOpRGJeHnZG38KB+/crLeNoaYlfRo1Eezc3LUZGiGpcLhfTQkegQVMXfL/gd0jEFZt37119iuKCElhY83UQof5gWK7SzdCMwoQYBhkjx9Jbx3AoIUZxLMjdB//rPAZmPOVfdWtT89IanfuX8MODy4rjbhY2+KXXOLSyL60GD3BphC/8BsDB3FIrz6G8PFHpZH+25qVv8ndTU/HzzSiEPX0KuYpOz280dXDAr6NHwdO2ZmtYEaItQeMC4ebphK+mbYUwu1DpnEQsxe2IWHQbYtw1izJwIGOxKCWbaw0BJTXEIIjlMnwWdQhhyQ8Vx4Z7tsfXHUfARMVQbpFMiiU3TuDoyweKY23sXbG95zi4WdoojnlZ28HL2k6jsZdXIpUg/OVz/BMXC7FMhp+DRyL82TP8EhWFG0nJFcq3dXXFs6wsFEtLv/l2adgA20aMgC3fuL/tEv3VNrA5Np4OxbJJm5H4NFXp3I0z94w+qSHsUFJD6r1iqQQf39iPS2lximMTm3TCF76qlzHIFhXhvUsHEZWZpDjW36M5vgscUaEjsDamc5fIZLic/BLH4h7h9IunKJRI4GRhibmtO2PIb7vxXMXcMv28vTG7kz+6NGyIjlu2olgqxXAfH6wJHghzE/qzJvrNo4kLNpz6DF/P+Akxlx4pjt88d9/ol1GQM+wm0JNXXolrFOjdj9R7sXkpuJb+XLE/p0V3LGzdX+UbY7wwC7Mi9uNluRFO01t0xucd+qucnE9T5AyD6NRk/BMXi5PPnyC7pLj0BANwJRzIJAzWXrykdI0Zj4eRrVthlr8/mjk6Ko6XSKV4L6ALQrp3pyHbpN6wsbPCqv0fYcviv3Dq99Lf9azUXDy7l4hm7b10HJ3uyFn2qWFzrSGgpIbUex0dvbCu8xiE3DiIj1v3xdyWPVWWi85MwtyIA8gVlyYQXA4HSzsOwJTmnbQSJ8MwiM3OwD9PY3H82SMkF+SXnZQDHDEHHCnAAQdCiBSn7Ph8TO7giykdOsDZykrpnnKGwfJ+/TC+fTsQUt+YmJrgow3voGEzV/yy/G8wDIMbZ+8ZdVJD2KGkhhiE4AatcSzofXjbVL7IZPnOvuVHOGlaiVSKX+5G4cjTh4jLVTFPhxzgFnHA+U8HPy9bW8z098eYtm1gaWqq8t5cDocSGlKvcTgcjPlgIDyauGDNu7/gxtl7mLRoqK7D0hk5OJCz6OzL5lpDYNz1VMSgVJXQAEATGwds6zEGXtb22Nd/ilYSGqB0vpg+Xk3gJbBTXYALoNxSUh3c3bBl2Fs4N3MGpvh1qDShIcSQBA7pgP8dX4yslBzkZuZXf4GBejOjMJvNmFFNDTEqnV28cGbIXJhytbsgZVsnV2wbOAJzTx/BhcT4Cue55hx0dmyIRd27w9/Dw6g7ShLj1ay9F74LW4LstDzYOdlUfwEh/0FJDdGp8LTL6OPSDTyO9ioNtZ3QAEB8Xg4+Oncc9zLTVJ7/smcfTGvbUctRGQapVAYTk/q1ajqpnJO7PZzc7XUdhs5QR2F2jPvZE51hGAZ/vDyM61kxWk1otI1hGPz95AGGHtytlNBYlxs6PtS7Jaa2obk56iIuIQPDP9qOn/++iqy8wuovIETPycFymQQj71NDNTVE62SMDD8924N/06/inUajdB2OxuSLRfji0jn8ExerOGbOM8GX3fpAJJVhZeS/8La1x9rewTppbpIzDDILC5EiLECKMB8pwnwUiESY2skPthb1Y/K+fadvITuvCL8eisTuozcwINAH4wf5oWVjWhiREGNESQ3RKpFMjI1PfkVUzl0AQFvbljqOSDNi0lPwUfhxJAjzFMdaOjhhU/+haOngjEtJL8A3McHWAcNhbab+lb8ZhkFucQlS8kuTlVRhgeLnN1tafoFivSsAaGgrwLa3R9SbhEYqk+P+0xTFvkQqw8lLD3Dy0gP4+TTE+EEd0dO/aYX5hwqLxbj7JBmBvk20HTIh1WJYjn5iqKaGEO0olBZh7aMfESssnfnXgsdHY6uGOo5KveQMg5/u3MD/bl6BtFzCMKV1B/xfYG/wTUpHMjW3d8TKHkHwcXRWewwMw2DzlevYdCmyxtf4N/TAltHD4Gil/TWu6sqEx8Ufa6bi8q1n2Bt2C7djy2aIvv0oCbcfJcHdWYBxwR0xrHdbWFuaAwBMTbhYtO4wQqb2xdsDqdmP6Be2K23TKt2EaEGOOA9fP/wBL4vK1i9qLWgOHsdwOnimFxYg5N9TuJz8UnHM1pyPtb2DMahJc6WyrpbWGNuyrUbi4HA4+LBHV7jZWGNpWLhScqXKqHatsXJQ/3q5vAKPy0XvTs3Ru1NzPHmRjn2nb+HM1UeQSGUAgJQMIb7/4wJ+PngFQ3u1xdhgP3i52YPL5WD9b+fxMiUHH7/TByY8w+3XReoX6ijMTv17FyP1TkpxOlY93IR0UZbS8Ta2LXQUkfr9m/Aci/49VbbcAYAu7g2xsd8QeFgLKpTXdB8amVwOOws+PO1sEZ9dce0oAOAA+KRvD8wJ6GQQQ8hbNHbBl/MG4f0JPXH43B38fe4OcoRFAICiEgkOnLmNg2dvo1sHb/C4XEhlchw4cxtJablYNX8orF7X5BBC6i9KaohGxRck4uvYH5AnqTiZVhtB/U9qRDIp1l6/hB33ohXHuBwOFvh3wwd+AVpdTwoAcotLcODOffx56w6S84SVlrM0NcX/hg9GUIumWoxOOxxtrTB7TDdMHd4F5649xt5T0XjyMgMAwDDAldvPlcpH3onH3BV7sf6TUXB3rpiAEqJN1PzEDiU1RGMe5D3B2kc/olhWUuGcFc8Cjep5fxqRTIoxR/7C/XJDtT2sbfB9v6Ho7K7d5/Y4PRO/R9/GP/cfoUQqrbKsh8AG294egVau6u/Po0/MTE0wpGcbDO7RGjGPkrE3LBoR0XFgVKxi/CwpE7OW/YlvQ0aibTN37QdLyGu0TAI7lNQQjbiedRsbn+yAlFH9AdvatkW9n5/GnGeC7g28FEnN4CbNsaZ3MGzNtTN6SCaXI/zpc/weHYNrLxMrnHewtMDkjr6w5fOx6twFAIBfA3dsHTMMTv9ZGNOQcTgc+LVqCL9WDbH94BXsOHxNZbnsvCJ8sGo/vnx3EIK6GuaoPEIMHSU1RO3yJQV4KHyKAMcOyBLl4FH+swpl2giaq7iy/lnUuQdi0lMwsnlrTPBpp5W+KdU1MTnbWIFjAvxfvz4Y0qIF9sXcAwCMaOODr4cMqJcdgtmSyeX4YU8E9p6KrrKcSCLFFz8cR0JKDmaMDDCIvkakfqHmJ3aM792NaJyNqTVmNBkHAPgr4R/VSY2BzE9jxuNh77DxWvnwq6qJicMBTEy5KIEMqZICfNKlB4a0KO2zVCgWY1Hv7pgX2NloP6SfvswAGAa9OzVDaqYQqZlC5BVUbBZ9Y/vBK0hMzUHo7AEwM6W3SaI9lNSwo5W/1i1btmDdunVITU2Fr68vfvjhB3Tp0kVl2QcPHmDp0qWIjo7Gy5cv8d1332HBggWs7kl0I6U4HUeTzyn2Ozv44mb2HVibWMHL0kOHkamXNhKFDRev4MerN1Q8OCDnMWB4gJRTOox5hI8P3uvcWVFkrG9b2Jgb98genyau8GmiPMtwcYmkNMHJEioSnZTMsp9PX4nFq4w8rF04HHY29Wf+HkKMmcaTmn379iEkJATbtm1DQEAANm7ciODgYDx+/BguLi4VyhcVFcHb2xtjx47FwoUL1XJPon0Mw2Bn/H5Fn5pWgub4pOVcrHy4CZY8C3DreX8abWIYBq4Ca+VjHAaMCcBwgfL9Aju4uWHNgIFKiZaxJzSVseCboklDRzRp6KjyvFQqQ3p2AQqLxZTUEK2hmhp2NP7JsmHDBsyZMwczZsxA69atsW3bNlhaWmLHjh0qy3fu3Bnr1q3DhAkTYF7Jm3Ft70m0LyrnLm7nPgAAcMHFrCbjweVwMdd7Ejraa2bSOUN16OFDfH3pAhgOAzmXgcyMgdwcYHhQSmjcra3x0/ARRtlnRhNMTHjwcLFFAxc7XYdCjAirxSxZJkSGQKNJjVgsRnR0NIKCgsoekMtFUFAQIiNrPoU723uKRCIIhUKljWiOSCbGzvgDiv1B7r3RyKoBAMDdwgX9XLrpKrR6aUybNjg/cyaGtG8Jxgwq/2otTEywfcRIONfTUU0nbsTi19M3kFdYeT8XQgipjkaTmszMTMhkMri6Krdlu7q6IjU1VWv3XL16NWxtbRWbp6dnnR6b1MyR5NPIeD17sK2pAOM8hymdN9bOqmzki8RIzMur9Pz6QYPQRkNNr2KpFMKSmicbkbEvsf7vi9h7MQaXHsQjPjUbIknVc+f4NWuAH45dQfCXP2PtgX+RlJnLMmpC6icGZXPV1GVTMQ2TUTGKeurQ0FCEhIQo9oVCISU2GpJakoF/ks8o9qc0GgUrEwsdRlS/SeVy/BwVhe+vRUIsk6ksszCwGwY3r/3szHKGQW5RMdLyC5CWX4D0/EKVPze0E2Dz2GEQ8Gs2/87NJ4n4499bFY672FmjoaMtGjqVbg2cbBX77vY28PV2x53nKfjrYgz2RdxBP9+mmNq/E9o3ocnwiPGgPjXsaDSpcXJyAo/HQ1pamtLxtLQ0uLm5ae2e5ubmlfbPIeq1M34/JK87B7e0aYpezgE6jqj+epyZicVnTuNeud91cx4PXRo2xKWXpYtmvtWyJeYH1Ow1joh7gYMx95H+OmlJzy+ApJrFLt9q64Ov3xoAfi2GNSdlqa5RSs8tQHpuAW49S65wjm9mApNyS0rIGQbnYuJwLiYOvt7umNrPH33aN9X6shOEaBslNexo9B3CzMwM/v7+CA8PVxyTy+UIDw9HYGCg3tyTqEdU9l3cyrkPAOCAg9neE6ipqQ4kMhm2XL+OEXv+VEpo/N09cOKdKXj39XDtdq6u+HbgwBq/xt29veBiY4XbSSlIzhNWmdBwACzq1wPrRw6qVUIDAGN7tEfIqF4Y19MX3Vs3RiMXe5iaVL0ae4lYioISscpzd56nYNEvxzHyq13YezEGxSJJreIhhBgPjTc/hYSEYNq0aejUqRO6dOmCjRs3orCwEDNmzAAATJ06FQ0aNMDq1asBlHYEfvjwoeLn5ORkxMTEwNraGs2aNavRPYn2lXYO3q/YD3brjcb1fG0nXXiUmYHFp0/jfnq64pg5j4dPevTA9A5+4HG5MOFy4Wplhe3DR4BvYlrje/O4XPzfwD7IKSrB8fuPKi1nZWaG/40ajL4tvOv0HDq38ETnFsrNuzK5HBl5hUjKzENSZh6Ss17/m5mHpKw8ZOcXVXvfxMw8rDnwL348EYmxPdtjQu8OcBLUz47RhFRG2zU1ERERWLduHaKjo5GSkoLDhw9j5MiRVV4jEonw1Vdf4Y8//kBqairc3d2xdOlSzJw5s85xq4vGk5rx48cjIyMDS5cuRWpqKjp06ICwsDBFR9+EhARwy1Upv3r1Cn5+for99evXY/369ejduzcuXLhQo3sS7UsuTlUsXCkwscYEr2HVXEHKk8hk2HbzJjZfv6ZUg9LJowHWDhyIJvb2imPuNjb4afgIuFpbq7qVSvklIhy5+xB7ou7ieVZ2peW87G3x4/gRaOaseu6WuuJxuXCzt4GbvQ06Na+Y7N58kog5mw7W6F5mpjzce5ECsUSGqUH+lNgQg6LtpKawsBC+vr6YOXMmRo8eXaNrxo0bh7S0NPz6669o1qwZUlJSIK+mKVtbOAyjas1awyYUCmFra4u8vDwIBAJdh2Mw8iWF2JtwFM1tGqOPCzUF1lRsRgYWnzmNB+VqZ/gmJvi0ew9M8/MDl0UT3pP0TPwZdQdH78aiSFJ1s023Jl74bsxQ2FloZ0HON0rEUkz69k88Ty1LtjgcoKGjLZq4OZRuro7wfv2zjQX1jyPapY3PjDeP0ePoBzCxqvvvuLRQhMvDt9QpVg6HU21NTVhYGCZMmIDnz5/DwcGhznFqilGMfiLaYWNqhTlNJ+o6jHpDLJPhxxs3sOXGdUjLfcvp3KC0dqaxnX0VV1d933OP4rAn6g5uJlTslCvgm6Nvc2/8cy9WcWxqFz98NqCXUmddbTkZFYum7o4I8msOb1dHNHFzQCMXe/DN6O2JGB+G4YBhUVPz5tr/zsemrgEzR48eRadOnfDtt9/i999/h5WVFYYPH46VK1fCwkL3I13pXYMQHWAYBpMPHkD0q1eKYxYmJvisZ0+849uhTrUzacIC7Lt1F/tv30NGQcU+Kq3dXDCpky/eatsSL7Ny8c+9WJhyuVg2pD/G+ululufR3dphdLd2Ont8QvTJm/lm2FwPoMK0JcuWLcPy5cvZhAYAeP78OS5fvgw+n4/Dhw8jMzMT77//PrKysrBz507W92eLkhpCdIDD4WBYSx9FUhPQsCHWDhgILzu7Wt2HYRhcf5mEPVF3cO5RHGT/aU025fEwuHVzTO7UAb4N3BQjpXJLSuBoZYkfxr4Ff88GanlOhBD9kZiYqNT8pK5pTeRyOTgcDv7880/Y2toCKF266O2338bWrVt1XltDSQ0hOvKOry8uvohHnyZNMLm9b51qZ8IfP8MHB45VOO5ha4MJHdvjbb+2cLSquBijvYUFDs6aCA9b6lNGiD5RV0dhgUCgkf4/7u7uaNCggSKhAYBWrVqBYRgkJSWhefPman/M2qCkhhAd4XI4+GXESFZz+fRq1hiOVpbIKixtburh3QiTOvmiT/MmVU5U19LVqc6PSQjRHHX1qdGU7t2748CBAygoKID16xGYT548AZfLRcOGup/Gg5IaQnSI7eSEZiYmmNG1IzIKCjHR3xdNHOvWuZgQYpwKCgoQFxen2I+Pj0dMTAwcHBzg5eWF0NBQJCcnY/fu3QCASZMmYeXKlZgxYwZWrFiBzMxMfPrpp5g5c6bOm54ASmoIqffmdOus6xAIIWqi7XlqoqKi0LdvX8X+m3USp02bhl27diElJQUJCQmK89bW1jh79iw+/PBDdOrUCY6Ojhg3bhxWrVpV55jViZIaQgghRE9ou/mpT58+qGq6ul27dlU45uPjg7Nnz9Y2NK2gpIYQQohBYhgGcbfjYWFjgYbN68dq7wzLmhpN96nRd7TkLSGEEINz/UQ0pjX/EO93+gyHvjuuOC6VSHUYFdE0SmoIIYQYHBtHG6Q8L13l/vLh65DJZBBm5+Pb6Zt1HFnVGAAMw2LT9RPQMWp+IoQQYnBaBTSHs6cjMhKzkJOWh6v/RGHP13+jIKdA16FVSQ4OOGqYUdhYUU2NmjFyYfWFCCGEaMSLB4l4eO0J8rML0GtMV8XxVeM3IO52PMQian4yZFRTo0aMvABMwRZwBKG6DoUQQoyStZ0lQnp9ifycQpiWWxRVLitdNFZSItZVaDWi75Pv6TuqqVEjpnAbIHuu6zAIIcRoOTVwxMfb5gEAJOKKtTISPa+peTNPDZvNmFFSoyaM9CVQuBOQ63d7LSGEGLreYwMRNKWXynNikUTL0RBtoqRGTZj81QAkAENJDSGE6Nr8TTPh2si5wnG5TA6ZVKaDiGqG1cin15sxo6RGDRhRBCA6/3qnULfBEEIIgZWtFT7b/aHK9dX0ubbmTZ8aNpsxo6SGJYYRgxF+XXaAmp8IIUQvtOvZCuM+HV7huESPkxrCDiU1bBX9Dsjiy/aZgirX0SCEEKI9074aj6YdGisd0+fOwlRTww4lNSwwsgwwBf+dnVIKQKSLcAghhPyHqZkpQv/4CGZ8U8UxsR4P66bRT+xQUsMCU7BBdR8aOfWrIYQQfdGotSdmr3lHsa/fNTXUUZgNSmrqiBHfAYr/ruQk9ashhBB9MmL+IPgP9AVAfWoMGSU1dcAwcjD5q6ooQEkNIYToEy6Xi092vA8bB2u9TmpKa1vY9KnR9TPQLVomoS5KjgCSO5WfpxFQhBCid5w8HLBg21yIS/Q5qaFlEtigpKaWGHkRmOKTgOVkcEyagxEuV1GIkhpCCNFHvd4OREkRDeYwVJTU1BKHawmOwy8AAKbk37ITph3B4Q8unVmYkhpCCNFbfEtzXYdQKeb1xuZ6Y0ZJDQuM5GbZjlkXcKymATwvgBHqLihCCCH1FjU/sUNJDRvisqSGY9a59F9+X5p8jxBCCNEBSmrqiJEXApIHr/e4gKmf4pyqtUYIIYSQalH7Eys0pLuuZEkA16H0Z5PW4HCtdRsPIdUQlpQgNZ/6exGi19gukWDkzU+U1NQRx7QlOM6XwHE6C45gma7DIaRSUrkcf0THYPKfB2FvwVc69ypPiN+u34JEJtNRdISQ8mhGYXao+YkFDocDmDQC0EjXoRBSAcMwuPjsBdacj8CzrGwsD+4HcxMTMAyD6MRk7L5+G2cfP8MPb78FUx5P1+ESQghrWqmp2bJlCxo3bgw+n4+AgADcuHGjyvIHDhyAj48P+Hw+2rVrh5MnTyqdnz59OjgcjtI2aNAgTT4FUgvUUVr3nmRkYua+w5hz4AieZWWjga0AI9r44NCdBxj9yx5M/u0ATj+Kw6BWzRHk00zX4RJCXqNVutnReE3Nvn37EBISgm3btiEgIAAbN25EcHAwHj9+DBcXlwrlr169iokTJ2L16tV46623sGfPHowcORK3bt1C27ZtFeUGDRqEnTt3KvbNzfV33gFjczXrJro7ddF1GEYpq7AIGy9dxf6Y+5CXSy6bOzoiaPNO5BQVK47ZW1rgy0F9dREmIaQybPvFGHlSo/Gamg0bNmDOnDmYMWMGWrdujW3btsHS0hI7duxQWf7777/HoEGD8Omnn6JVq1ZYuXIlOnbsiM2bNyuVMzc3h5ubm2Kzt7fX9FMhNfCs4AXOpV3UdRhGRySVYvu1mwj6aSf23r6nlNCAAS4+iVdKaADgy+C+cLCy1HKkxo1hGMjlVJNJiKZoNKkRi8WIjo5GUFBQ2QNyuQgKCkJkZKTKayIjI5XKA0BwcHCF8hcuXICLiwtatmyJ9957D1lZWZXGIRKJIBQKlTaiGSdSziK5OEXXYRgNhmFwKvYJBv38G9b9exkFInGFMhw58N/vbv1bNMWQNi20E6QBYhgGUbfia5ygyKRynD/3AO/O/hUXzj/UcHSkPqOOwuxoNKnJzMyETCaDq6ur0nFXV1ekpqaqvCY1NbXa8oMGDcLu3bsRHh6OtWvX4uLFixg8eDBklYzgWL16NWxtbRWbp6cny2dGVMkoycT1rGjkSwsglOTrOhyDd/dVKib+sR8fHTmBpNxKEnUVc14I+OZYNqQfzadUR4+epGDB4r/wyef7cfb8g+ovAHD0SDS+WfkPnsWlY99fkdTvjFSOUcNmxOrlkO4JEyZg+PDhaNeuHUaOHInjx4/j5s2buHDhgsryoaGhyMvLU2yJiYnaDdhInEoNhxxyAMCrYtVJK1EPhmGQXlAIL3s72PH5lZZTVUsTOrA3XG1oXqW6uhEVjzv3St9DftkVgZIarPg8ILgdLC3NAADP4tJx88ZzjcZIiLHSaFLj5OQEHo+HtLQ0peNpaWlwc3NTeY2bm1utygOAt7c3nJycEBcXp/K8ubk5BAKB0kbUq0BaiPPplxX7r0ooqdEkDoeDoBZN8e1bwfhm6ADVhVR8a+vRtBFGtW+t8fh0iWEY/HH0Bv44drP6wnUwbnRnODmWJoUZmfk4eCSq2musbfgYNqKjYn/fHtXN74TQ6Cd2NJrUmJmZwd/fH+Hh4Ypjcrkc4eHhCAwMVHlNYGCgUnkAOHv2bKXlASApKQlZWVlwd3dXT+Ck1sLTIiCSixT7yVRToxWX41/i4yMnVZ77by2NpZkpVg4NMuhmJ4lUhlXbTmPLX5ew9a8IRESp/qLDBp9vilnTein29+y/hpzcwmqvGz22C0xNS+cDuhOTgNiHyWqPjRgIanqqM403P4WEhODnn3/Gb7/9htjYWLz33nsoLCzEjBkzAABTp05FaGioovzHH3+MsLAw/O9//8OjR4+wfPlyREVFYf78+QCAgoICfPrpp7h27RpevHiB8PBwjBgxAs2aNUNwcLCmnw5RQSKXICz1vNKxV9RZWONuJCThvYNHFbMBO1pagPPmTY0BbEzNsH3CSPBNSmduWBzUEx62hl1LyQGQkV3an4thgOWbT+Lpy3S1P87A/m3Q1Lt0SoqiIjF+++NKtdc4OlojaGA7xT7V1hCifhpPasaPH4/169dj6dKl6NChA2JiYhAWFqboDJyQkICUlLIPwG7dumHPnj3Yvn07fH19cfDgQRw5ckQxRw2Px8Pdu3cxfPhwtGjRArNmzYK/vz8uXbpEc9XoyNXMm8iV5Ckdoz41mhWTnIK5B46gRCoFAFiZmSJbWIzXXZrgYmmJ/bMmonfzJvBt6I4ujRpifMf2OoxYO0xMeFj18TB4uZdO8VAskuDTdUeQVYOalNrg8bh4b3Yfxf7RkzF4mVj5CMw3xk8MwJuKsiuXnyDhZaZa4yL1HzU/scNhjLAbvlAohK2tLfLy8qh/DUsMw2Dx3RVIKn6ldJwDDnZ1+QFmXDMdRWa4HqalY8qegxCWlDb3uVhbQSqWIbeoBAwAO2s+wuZOU8xB8/PVmwj2aQ4vBzvdBa1liSk5mPXlHuQXlgAA2jRzx5Yvx8HcrLTWSlhQAoF15R2sa2rxF/txIyoeANA9sDm+Xja62mtWLP0bly4+BgAMGuKLTz4byjoOolna+Mx48xie25aBa1H33015cQkS311htJ9v9XL0E9Efd/IeVEhoAIABg5TiNBVXEDYScnIx/a9DioTGycoSf0wai23jR8CMx8PYDm1x+t3pSpPqzejqb1QJDQB4utvjmwXDwOOVvsU9iEvBN9tPK4ZSb917Cc8S2deSvDu7L7jc0m/GVyKfIuZuQrXXjJ9Y1j/w3Jl7yMyg6Q9IeRw1bMaLkhrCyolXZyo9R01Q6ucmsEHXRg0BAPYWfPw2cQyaONrDr6EHjsyZjFVvBcHewkLpGhOucf6Zd2rrhUXT+yn2z1x5hN+OXEdefjFORTzEwdO3WT+Gd2NnDC7XT+bHn/9VTMiXm1uE9IyK8wf5tPKAX8fSRXClUjn+PlC2Fl5BfgnrmAgxZsb5bkfUIr4wAfeFjyo9TzMLq58Zj4cNI4Zgckdf7JwwGi2cnRTnmjo7GvTIproYFeSLcYP8FPs/7b+CZZtPQCyR4tSlhxAWsE8iZk7pCb65KQDg8dNU/HsxFgzD4H+bTiPmjuqam/GTymprjh+9jfz8YkgkMnwRuh9ikZR1TKQeo8n3WKGkhtTZxYyr6OfSE//XaiE8+GXzCLW3bQOAhnVrigmXi+XB/dDGzbX6wgQfvtMHXX0bK/av330JABCJpTh24R7r+zs6WmP822ULuG7feRHHT93BpatP8PBRxaZZAPDv1ATNmpf+/ysuFuPYkVvYsukM7t9LQkpKDuuYSD1GSQ0rlNSQOpvaaBzmeE9BG4EPMsVlIz8WtpiHRS3eR5GsSIfREVJKLJFieN92iv415f19JgYyuZz1Y4x/uwsc7K0AAGnpQvxv02kAwMPHqmsrORyOUm3NH7uv4PjR0uawpMRs1vEQYqxMdB0Aqb+4nNIPiVxJHsTy0qnibU0F4PP46OTQAa0FLXUZHjFyDMPgl4NX8duR65BVsvBkSoYQV249R69Ozer0GE+fpeFpXBpKSiRo3MgJ2TnKQ8efPU9HSYkEfH5p85REIsPvv12CRCwDl8cFl8eBXMZALC5rckpKopoao8ZwSjc21xsxSmoIa2klGYqfXczL+nhYmlioKk6IVnA4HMx+uxs8XGyxcfcFFBSJVJbbH3a7zkmNm4stNmw6jdhKamRkMjmexKWhfdvSzt2mpjz4dWyM0E/3QipVXUOUnEQ1NcaM7UrbxjdJizJqfiKspYnKkhpXvrMOIyFEGYfDwdDebbFn3XT06Oitskz0g4Q6D++2seFj3Tfj0bZ1g0rLPHykvByCX8fGWBw6rNLySZTUEFJnlNQQ1tJLyj4QXM0pqSH6x9nBGt9+MhLLPxiictK9AyyGd1tbmePbVWPh285T5fnYRxVrcfoFtcG89/urLE81NUZOyx2FIyIiMGzYMHh4eIDD4eDIkSNVlr9w4QI4HE6FLTVVPwaGUFJDWCtfU+NCNTVET3E4HAT3aIW/1k1H3y7Nlc6FsRzebWlpjjUr30bHDo0qnKtsBNTY8QEYM7ZLheNZmQUoLhLXORZSz73pU8Nmq4XCwkL4+vpiy5Yttbru8ePHSElJUWwuLi61ul5TKKkhrKWXUPMTqT8c7KzwzcLhWPXxW7AXlPb7Usfwbgu+GVavGIPOHRsrHc/IzEdGpupZg+e93x+9+7aqcDw5mToLE+0YPHgwVq1ahVGjRtXqOhcXF7i5uSk2rp5M8qkfUZB6TalPDTU/EQ2QyuSIeZJcfcFa6N+1Jfasm44B3XwAAAdPsx/ebW5uilXLx6Brl6ZKx2Mrqa3hcjn47PNh8O3gpXScmqCMF4dhvwGla0mV30Qi1R3l66pDhw5wd3fHgAEDcOVK9avUawslNYSVYlkJhJLSb6HmXDPYmhrfAmqkzOaTV1Eslqj1nskZeZi3Zj/e/fYAHjwva7c/d/MJ63vbCSzx1YdDsXbRCEikMlyOfsb6nuZmJlj55Sj06FbWxFVZExQAmJmZYMWqt9G4SdkXAkpqjJia+tR4enrC1tZWsa1evVot4bm7u2Pbtm34+++/8ffff8PT0xN9+vTBrVu31HJ/tiipIayk/2c4N03Tb9yuPnqBCev/xMNE9S1m+r89/+Ju3CvIZHIs3X4SRSViPE3MwPKfT1U6TLu2enVqhj3rpqNATX1ZTE15WP75CPTpWTpXU1VJDQBY2/Cx+tvxcHa2AUAjoIyamvrUJCYmIi8vT7GFhoaqJbyWLVti3rx58Pf3R7du3bBjxw5069YN3333nVruzxYlNYSVdFG5kU9G1p9GIpdBLJPpOgy9sXDHMaTnFeBFeg6mfLcXP5+5DqmM/Wy9i9/pDxtLcwBAQlouvtt7EVv/vgyRRIaIGPY1K28IrPkY2ruN2u5nYsLDF0uGI6hvazx+mgqptOrfFWcXAVavmwBraz6SaQI+wpJAIFDazM3NNfZYXbp0QVxcnMbuXxuU1BBWlCfeM66k5n5mGi4mxus6DL2QlJmL83fjkJFXOqOuVC7H5pNXMWvzfiRl5rK6t5ujDZZMDVLsH7l4D5fvlL7uZ2+wb4LSJBMeF6GfDEWfnj6If1H9XDiNmzjjq2/eRnpanhaiI3qpHq79FBMTA3d3d+0/sAqU1BBW0ox45NPNlCQcevpAa493NyUVcj2dLtRRYIU+bStObhcTn4Kx6/7AkesPwLCIvVu7xujYsmGF49fuv0B+Ue2HYgvzi3Hs3F0kvNJ8Mw+Px8VnIUNgZ2dZo/Ltfb3w3ocDUFio3o6dpJ7QclJTUFCAmJgYxMTEAADi4+MRExODhITSFeZDQ0MxdepURfmNGzfin3/+QVxcHO7fv48FCxbg/Pnz+OCDD+r6jNWKlkkgrMgYGUw5JpAwUqObo+ZGajIiEuORJyqBrXnFCd3UbcfNW/ioe1d4Ozpo/LFq65/rD3Dh/nOV54pEEiz76wwuPniOpeOCYG9d8+Uzoh8lYun2U0jPKVB5XiqT4+KtZ3irR82bjfYdi8LW3yMgk8kxZXQA5k3uWeNr64rL5cDZyabG5Xv19mGVBBJSU1FRUejbt69iPyQkBAAwbdo07Nq1CykpKYoEBwDEYjEWLVqE5ORkWFpaon379jh37pzSPXSJwxjhX45QKIStrS3y8vIgENBoHbbkjBw54jxYm1jCnKe5dlt9ImcYdNy9GbmiEnzTcyAmtfLV6OPlFpeg25btWBncH2Paqa/fh7rkFhZDLJVBLJGW/vt6k8jK/SyVQWDJR5fmqmferczNhwn4v20nkJNfrPJ89/ZNsHFhzefYiLz1HJ9+fQgA4OFqi31bZlMHd1IlbXxmvHkMz/UrwbWo+5ckeXEJEj/50mg/36imhrDG5XDhaG6v6zC0Ki4nC7mi0maPQ08eaDypORb7CGKZDDGvUvQyqbGz0tzipZ1be+H35e9gyZZjuP+84lTs1x68hLCwBAKrmn0QdG7fCAJrPoQFJXiVlofHz9Lg08xN3WETUje0Sjcr1KeGkDq4kZqk+DkqLRkJwlyNPt7Bu6V9d24nq14N2tC5OtjgpyXjMKZvxeRRJpPjwq2aj7wwMeGhd9cWiv3wK4/UEiMhpHZ+/PFHtd+TkhpC6uBmuaQGAA4/faixx3qYlo4HaekAgCeZWSgQGee6QGamJlgytT+WzQqGmQlP6VxtJ+Lr372l4ufwq48hlxtdKzzRU+qaUbg++Pnnn9V+T0pqCKmDqFTlKfsPPWU3uqcqf98rG2ElZxjc0+JquKl5qtcs0iVbKz5kYpnSKI8bDxOQW1Bc4+SkQxtP2NuWjkZKz8zHgydVT45HiNbUwyHd+oSSGkJqKblAiOQCodKxl8Jc3EpX/wejSCrFPw+Um0duv9J8UpNdUIRV/5zHwRvsFnnUBN8WDWBqwgNHDsUbuEwmx9rfziExrWaT1pnwuOgbWL4J6jGA0qHex8P17zkTYogePHiAwMBAzJ49G9999x3Onj2LV6/YvY9SUkNILd1MSVJ5vK5NUHFZWZWeC497jtwS5XlYbidrrlZBJJHi14s3MXj9ThyPeYQpPTqq9f6345JRWMKu+UxgxUcf/2bgAODIAd9mHgCAczefQlhQ8zlr+nf3Ufz879XHyMsvxsKVB3HvkXoXziSEqNayZUts374dQUFByMrKwtatW9G7d284OzujZ8+6TbVAo58IqaUbqaqTmuPPHuHLwL4w59Xuz+rbiEv435DBsFExjfmbDsLlxbxKBcMwah2GzDAMTt97gg1hl5GcU1oL9eGAbrBlMbS0vOepWQi/HYcLd59hSn9/DOrUsvqLqvBWr7Y4c+0xOACeJ2SCxwAylC5+2a65R5XXymRy8HhctPNpAGcHa2RkFyArtxAzP9mNtMx8NHQzrpF8RL9wwK5fTH0a+8Tj8dCuXTu0a9dO6XhJSQliY2PrdE9Kagippf92En4jV1SCfxOeY1CTFirPV+Z+ehr+iLmD9wK6KB1/JczHpfgXFcrnFBfjZW4eGtvbKY4lZOfCy8GuQtma4nA48PXyUKzVZG9lgSnd/ep8v/IYhsHnO0/hcVLp7NPnbj+pU1JTIpZgyaZjKBZJ4CAom523sLi05ocDICG1+uanPw7fwLFzd9HAzQ5cbtlHQFpmaf8hsVha69gIURsjGtL97rvvqjzO5/Ph51e39x9qfiKkFnJKivE0p/Lmoto2QZVIJEgrKMSO6GgUSSRK547cf1hpn7+YV2VDuwtEInz375VaPa4q7nY2eK9/VwQ09cTs3p1hZW7G+p5JmbmYuOZPPE0uW/foyoMXKBZJqrhKNb6ZKT4Y1xOx8Wk4f/NphfMcQLHwZVUmj+oCZ0cbRN9LUCQy5ZWIax8bIaR6WVlZ+OSTT/D+++/j7t27mDdvHgAgKSkJBQWqZw2vLUpqCKmFm6lJMOVyMbd9ZwR6eCmOf9d3CEIDeuNpTiZyS1TPfKtKorC0qSe7uBj77pZ1UJUzDA7cq3xdqfLz1TxKy8CZR3FIz2f/pjC2SzsseasPJnRVz2SCDZ3s8M2MwfBysVMcK5FIcflB3RYCbe7ljM+m9a/0fF4N+tSY8LhY+vEQWFmqTtpEVFNDdMmARz/Nnj0b27dvx7Vr19CzZ09cv34dHTp0QKNGjeDo6IiFCxeyfgxKagipBVMuD2fGzsTnXfvAzdJacZzL4WKebxecHTsTZjxeFXdQlpCbq/j556ibEElLP1BvJiYjMTcPfh7uaOnspCjjblO6flD5mpr7KemQyuU4GKOexTVbuDmBb6q+lmlvN0fs/nQiercrW/DyXEzFmpaaGtqzDUb0bqvyXE07Cru72GLR3AEqz5WIKKkhOmTASU1ERAT+/vtv3Lp1C+vXr8eoUaNgZ2eHw4cPY/Xq1di1axd+//13Vo9BSQ0htdDXyxtNbEs7ktqYlTV15ItLV1TmcbmwNK15s01CXp7i57SCQhx6WNp8VSKR4NiMd3BgygR4CMoWQlwV3B9bRr6FfJEIRa+bSR6mlE7Mt//WPcjk8jo+M82ysTDHhrnD8d5bgeBwgEv34lHCokZk0ZR+aNnIpcLxvMKaj34a2LMVgnu1rnBcTM1PhGhETk6OolPwtGnTkJaWhrVr12L48OEICQnBhg0bsHXrVlaPQUkNIXWkKqmprfI1NQDw042bkMrl6N20CVq5lK56nlNc9kFtb2mB4JbNcWLmFDCvv5I9TC1NalKE+bjwtG7NOtrA5XIwd3BXfP/uCPB4XETGvqjzvczNTLD6w7dg/Z8+NMKCypv+XqXl4rtfw/HF+qPY/NsFAEDInP5wd7FVKkc1NUSXDH1GYS63NO0wMzODpaUlnJ2dFed69+6Np0/rXosLUFJDSJ3ZmJXVyNQ1qUksV1PzZv/YI+XJ9nKLyz6o7V4PsbYwNYWVmRmKJRI8y8xWnP8r+m6d4tCmnm298cfiiYhPy66+cBUauNhh2dxBSseqqqkpKpbg4Knb+PfaE9y89xIAYGVpjqUfDwGv3Cgo6lNDdMqAm58AYM+ePYiJiYFUWvHvzMrKCjk5NZtAszJaSWq2bNmCxo0bg8/nIyAgADdu3Kiy/IEDB+Dj4wM+n4927drh5MmTSucZhsHSpUvh7u4OCwsLBAUFsc7uCKkt9dTU5FU49uP1G5CXW3Ihr9zke/YWyqthP0rLUCp7+dkLJObk1ikWbWrkYo9pQZ1Y36dXx6aYOrSzYr+qPjWO9laKn7NyChU/t/NpgGljAxX7JXUYmUW0LyU+HfH3EhB3+wVExca5Hlp907NnTyxbtgz+/v6wsbFBUVERli1bhp9++gk3b95EUVER68fQeFKzb98+hISEYNmyZbh16xZ8fX0RHByM9PR0leWvXr2KiRMnYtasWbh9+zZGjhyJkSNH4v79+4oy3377LTZt2oRt27bh+vXrsLKyQnBwMEpKat6eTghbbJMaOcNUqKkBgGfZ2ThdLkm/8sFcXJs/F6dmTYW1mXJ/nTf9ad5gAOyNrh/T/PO46nn7mfd2d3T0aQig6tFPtjYWMDEpfcycvCJIpTLFualjuqJdy9JJ+8RiqcbW8SLq8/WkTZjnvwTvB3yOpCcGtHq9AdfUXLx4EXl5eXj06BF27tyJRYsWISUlBaGhoQgICEDz5s1ZP4bGk5oNGzZgzpw5mDFjBlq3bo1t27bB0tISO3bsUFn++++/x6BBg/Dpp5+iVatWWLlyJTp27IjNmzcDKK2l2bhxI7744guMGDEC7du3x+7du/Hq1SscOXJE00+HEAXlpKb23xTTCwoglslUnvvx+g3FB6sZjwcnKys0d3KsMIvwg9SKXw7+vvMAYhVVu4bKhMfFqveHwtHWCkUlYqVkpTwulwNHu7Lamuy8sm+FJjwuvnw9zFsmZyCV6meHa1KGZ1L28SWr5P95fWTofWoAoHnz5pgwYQK+/fZbnDt3DtnZ2Xj27Bn++usvfPbZZ6zurdGkRiwWIzo6GkFBQWUPyOUiKCgIkZGRKq+JjIxUKg8AwcHBivLx8fFITU1VKmNra4uAgIBK7ykSiSAUCpU2Qthi26cmQUUtzRsP0tNxIb76Tr//rakBgJyiYpyOjat1PPWZo50Vvv5gKHhcTpX9ahzty4bhZ2Yrz+vj4WqHRXNK31eoX43+45mUTZ1gSEmNsWrSpAnGjh2Lb775htV9NJrUZGZmQiaTwdXVVem4q6srUlNVrzScmppaZfk3/9bmnqtXr4atra1i8/T0rNPzIaS88jU1wrokNa/705ib8MAtVwPzaY8e6NygAbbduFllM4hIKsXTDNWzG/8VfafW8dR3fj4N8d7YHlX2q3Eqn9TklCU1b5KYgb1aY2CvVjpJahJeZODF83S8jM+ATEY1RaqU/3vgmVZMauRyOST1PSF9s0wCm82IGcXop9DQUOTl5Sm2xMREXYdEDADbPjVxWVlo6eSEI5Mno5GdneJ4v6be2DthPLYMewvSKuadeZKeCalcDtNyfVMcLC3QwtkR0Ymv8Dgts9JrDdU7QzrBxcFG6VhmToFiTSun/3QWzsjKR+i3/+BFUllyGDI7CGZmNZ9AUV0WzvsNcyf/hDmTtkFUQp2VVfl740nsX38MRfnFSjU1UokMN07dxoJeyyGp7x29DbhPjTZodEFLJycn8Hg8pKWlKR1PS0uDm5ubymvc3NyqLP/m37S0NLi7uyuV6dChg8p7mpubw1zFCsiEsME2qeng7o70/AI0d3REUwcHxL8eyvgsOxstnJzgZGVV5fUPU9Ph7+mBr4YGYervB5FVWAQnayscnfsO7qekKQ0FNxYcDgdWFsqdqVMzhJj56e8Y1Ls1xJKyb/EnLzzA1j8iIJZIsWLhUMVxayvdvFeUr4VQ5wrshqRV1+ZY2Hs59q07ChOzso+v9bN/QnpCJnz7tIaljUUVd9B/bPvF1Ic+NZqk0ZoaMzMz+Pv7Izw8XHFMLpcjPDwcgYGBKq8JDAxUKg8AZ8+eVZRv0qQJ3NzclMoIhUJcv3690nsSoglWpmZ489FTIBHXesTMoBbNcTE+HhsuX0FTBwfF8bisms3fEtjYC39MG4dmzo44OGsibn32AY7NmwIOh4N2Hm4IaEzNrADQprk7HO2t8Oc/N3H8fNkoyodPU1BULIa3pxPM1LgsRF0pJTVcSmpUaRXQDA5udsjPKUROWlmftPSE0lrJzgN9dRUa0RMab34KCQnBzz//jN9++w2xsbF47733UFhYiBkzZgAApk6ditDQUEX5jz/+GGFhYfjf//6HR48eYfny5YiKisL8+fMBlH6DWbBgAVatWoWjR4/i3r17mDp1Kjw8PDBy5EhNPx1CFLgcDo6Nnop/x8/G9cnv1ekejezs8OO1G8goLJs35Xl2zZIaLwc7RV8cD1sBrMzYr6ptiDgcDkYFd6j0fIsmFZdb0AVGXpbUcKmmRiUul4tuw/0rPd852ACSGmp+YkXjX0/Gjx+PjIwMLF26FKmpqejQoQPCwsIUHX0TEhIU0yYDQLdu3bBnzx588cUX+Pzzz9G8eXMcOXIEbduWLWC3ePFiFBYWYu7cucjNzUWPHj0QFhYGPp+v6adDiJK2Tq7VF6pCY3t73E1Nw5EHD/Gm2udZDZMaUnMDevhgy+6LKCiq2EzYogm7/4fqIqeamhrpNqIzjm8Pr3Dc0cMejdsaQO0k22HZlNRo3vz58xU1Lf914cKFCsfGjh2LsWPHVno/DoeDr776Cl999ZW6QiREJ7xedxCWM1BKauQMQ9/W1ciCb4bBfdrgwMlbFc619NaPpKZ8TQ39r6+cb+/WsLazREGu8uyznQb6Ul8kYhyjnwjRV43t7QC8zmcYwIzLw8T27SCpZFI+UncjVfS34HI5aNrISQfRVFS+SxZ9OFfO1MwEAUM6VjhuEE1PADU/sURJDSE61NjeXvEzhwEkUhnCnz5Hvqhua0mRyjVu6IiO/2me8PKwhwVfP/oile8ozKXmpyp1H6G8bhiXx0XH/m0rKV3PUFLDCiU1hOhQo9c1NeUl5uVhzqEjKBLX8/k29NCogR2U9vWlPw0AyJWanyipqYr/wPYw45sq9lt3bQ5ru6qnQCDGgZIaQnTI3sICtvyK86LcS03DR8eOVzn5Hqm9Xl2aKa3/pC/9aYCymhqqpamehRUfnco1J3YyoKHcxrD2kyZRUkOIjpWfTbi8C8/jsfTsOVoxWo1MTHgYFtROsd9cT4ZzA+VqaiinqZHyTVCdgtvrMBKiTyipIUTHyver+a/9d+9jc+Q1LUZj+EYM8AXvdW2IrueoEYtKZzhW6k9TrulJRquFVypgaEdweVzYuQjQrENjXYdD9ITup9EkxMip6ldjZWaGf6ZOBo/DhUgmpSHeauTiaIPunZoi7mUGbKx0O7fV10v/RvsOjTBkZNloHg6XA7FIilNHb+Ph/SSErhilwwj1l8DBGr59WsPRzU5prrN6j21nXyOv2KWkhhAde9P85GlrCy6Hg5e5uSgUixER/wJTO/rpNjgDNSq4A/45q/uVzFv4eOCnTWdx4M9IxTGpRIapb/+A7MwCzP0wSIfR6b/uIzrD2tZS12GoFa39xI4BpbeE1E+N7e1hYWqCbaOGY06Xsn4Cv96Mpo7CGtKpXSMEdffRdRgI6NYMAJCdVaA4xjBAdmbpfmDPFjqJq77oPrwT/Ae0q75gfUPDueuMkhpCdKyRvR3WDApGS2dnjGrTGk6Wpd88k4VCnHz0WMfRGSYul4M+XXWfMDRt4QZHJxuV57waO6GBp6OWI6pfHD3sYesk0HUYRI9QUkOIjtlbWGCoT0sAgLmJCaZ3Kutf8fONKBr9pCH6MBcMh8NBl9e1Nf9FtTRGiibfY4WSGkI0SM4wOPPkaa0Sk0m+7RUrbsdmZODSi5eaCo/ogc6BqpOabr1aajkSog9onhp2qKMwIRrE5XBwIykZ++/dx9fBA+BqbV3tNQI+HxN92+GXm9FwsrSEsKRELbHkFpXg34fPYG5qAgtTE1iYmcLSzBQW/91MTfSiFsNYdOzcBCYmXEjLDd+2d7BCy9YNdBgVIfUTJTWEaNg7fr4I+mUnBu/4DUv798WI1q2qTRqm+3dEY3t7jGrTGuYm6vkztbPkw0VgjdD9YcgqKKqybN9W3tgw+S2YmfDU8tikcpZW5mjn1wi3b8YrjnXt0YJmFjZWNKSbFWp+IkTDGtvbo0fjRhCKRPjkZBjeO3IUmYWFVV7jZmODCb7t1ZbQvNG9RSMcXjAFvVo2qbSMr5c71k8aSgmNFgV0a660T/1pjBc1P7FDSQ0hWvCOXwfFz+finmHwzt9w4tFjnXQCdrS2xNbpI7Dkrd4w5VVMXO4kpODqE+rHo00B3cv61ZjzTeHXqfKkkxBSOUpqCKnEo6wM5JQUq+Vefb2bwENQNnQ3p7gEHx87gQ+PHkdWUdVNQZrA4XAwpUdH/PXBBHg7OyidMzfhIaCZp9ZjMmYNPB3RwLP0/4N/F2+Yl1uBmhgZGv3ECiU1hFTC0cISw/b/geiUZNb34nG5mOhbcSXhsCdPMXjnbzj95Cnrx6iLVh4u2PfhJLzdpS0AoGNjDxxeMAVW5ma1vleJWKru8IxKl9dNUDTqychRUsMKJTWEVMLZ0gre9g4Yf2Q/fr7Nfr6Yce3bqmzuyS4qxgf/HMPC4yeQU6yemqHasDQzxYrRA7Bh8lAUlIjRyKnyBTYrc/xmLL4/dkkD0RmPgG7NwOVWPm8NIaR6lNQQUoWRLVpBKpfj66sXMefkP8hjMbza0dISQ1pW3gH0WOxjDNm5G5fiX9T5MdgIbtcC22aMqlXyll8sQujuU/i/P8Jgb21Ya/BoWzu/RujUtSns7K10HQrRIW13FI6IiMCwYcPg4eEBDoeDI0eO1PjaK1euwMTEBB06dKjdg2oQJTWEVGFgk2bgvx6BdO7FMwzd/zvupKXW+X7v+FVsggKABgIBlvbvi5/HjERgI686358tV1vrGs9Rc/t5MsZ9+wdORj8CACRm5uL3f29pMjyDZmrKw9wPB+g6DKJrWm5+KiwshK+vL7Zs2VKr63JzczF16lT079+/dg+oYTRPDSFVsDYzw4AmTXHsaekaTEn5Qrx96C983q03prf3q/UkdR3c3dHG1QUP0tKVjqfk56ORvR3aurqqLXZNkcrk+PnMdWw/fR3ycrU6R288BN/MBGO6tYOlOXV0rQuvxk66DoHompbnqRk8eDAGDx5c64d59913MWnSJPB4vFrV7mga1dQQUo2RLVop7Uvkcqy4/C/eP30MQpGoVvficDhKw7u97GwBlC6nsODYSTzPzmYdryYlZeVh5qb92BZ2TSmheaNELMXt5+w7VhNC2BEKhUqbqJbvVVXZuXMnnj9/jmXLlqntnupCSQ0h1ejl2Rj2fH6F46eePcWw/X/gfkZare73lk9L2PLN4ePshIOTJ6Ghbekqw/kiEeYd+kdtyyKo2/GbsRi39g/ceZFSaZl3evuhe6vG2guKEAOjrj41np6esLW1VWyrV69WS3xPnz7FkiVL8Mcff8BEzZODqgMlNYRUw5THw9BmqofZvhTmYvTBv/D7/Zgad7C1MDXF223bYmYnfzhYWuCnUSNgaVraXBOfk4MFx09CJpdXcxftKd8ZuFAkrrLsHxdvY1d4lJYiI8QAqalPTWJiIvLy8hRbaGgo69BkMhkmTZqEFStWoEUL/Zz1Wv/SLEL00MgWrfDH/Tsqz4nlMmyNvg4+zwRjW7Wt0f2m+XeEo6UFAKClszP+N3Qw3jtyFAAQEf8C6yIuYUmf3uoJnoWUbCG+PXwRWcJCBLTwhLmpCfhmpuCbmpRuZqalC2Sambw+ZwK+qSmKRGJY1mGuG0KIeggEAggEArXeMz8/H1FRUbh9+zbmz58PAJDL5WAYBiYmJjhz5gz69eun1sesLUpqCKkBfzcPNLQRIClfqHTczcoax8dNgZNl7YYzl59dGAAGNG+GhT264bvLVwEAv9yMRksnZ4xq25pd4Cy5Owjw3axhOo2BEGPCdv0mTa79JBAIcO/ePaVjW7duxfnz53Hw4EE0aaL75T2o+YmQGuBwOBjxnw7DAJBaWIAt0dfV8hjvdw3A4HLz2Hx+5ixuv3qllnsTQuoJLQ/pLigoQExMDGJiYgAA8fHxiImJQUJCAgAgNDQUU6dOBQBwuVy0bdtWaXNxcQGfz0fbtm1hZaX7OZYoqSGkhkY09wEAOFlYYlmPvorjO+/ewtn4ONb353A4WDsoGK1dXAAAplwusou0P8MwIcR4REVFwc/PD35+fgCAkJAQ+Pn5YenSpQCAlJQURYJTH3AYXSwTrGNCoRC2trbIy8tTe5sjMWyD9+7GuFZtMb29H+afOY4TcU8AAHbmfJyaMBXu1jbV3KF6r4RCfHj0OL4eOAA+Ls6s70cIYUcbnxlvHqPV+9+AZ15xtGVNyUQliN36udF+vlFNDSG18J5/F0xs0x4cDger+wxEQ5vSN41cUQk+PqueUUseAgEOTp5ICQ0hRoijhs2YUVJDSC0Mb+6jWDZBYG6OHwa+BRNu6Z/RjVdJ+CHqmloep7YzFRNlOcIifLzyAJ4nZuo6FEKIFlFSQwgLfm7u+CSgu2J/U9Q1XEtO1GFERC5nsHLzKdy8l4BZoX/izOVYXYdESM1puaOwodFoUpOdnY3JkydDIBDAzs4Os2bNQkFBQZXXlJSU4IMPPoCjoyOsra0xZswYpKUpz9jK4XAqbHv37tXkUyGkUnP9OqOnZyMAr5c7OHsS2cVFOo7KeP159CauxbwAAIglUjja6X5EBiE1pe1Vug2NRpOayZMn48GDBzh79iyOHz+OiIgIzJ07t8prFi5ciGPHjuHAgQO4ePEiXr16hdGjR1cot3PnTqSkpCi2kSNHauhZEFI1LoeDDUGD4WRROldNamEBFp8/XeMZhon63HucjO17Lyv2Z4wJhH9b3a16TkitUU0NKxpLamJjYxEWFoZffvkFAQEB6NGjB3744Qfs3bsXryqZeyMvLw+//vorNmzYgH79+sHf3x87d+7E1atXce2acl8FOzs7uLm5KTa+irV5CNEWZ0srfDegbKXbcy+eY9fd2zqMyPgIC4qxdOMJyOSl7+od23hixttddRwVIUSbNJbUREZGws7ODp06dVIcCwoKApfLxfXrqicri46OhkQiQVBQkOKYj48PvLy8EBkZqVT2gw8+gJOTE7p06YIdO3ZU+a1YJBJVWLGUEHXr6dkY73XsothffTWi1otdkpq7+ygZwoLSxT8ZhsGqraeRlpUPALCzscCyD4eAx6Vug6QeolqaOtPYX3xqaipcXk8i9oaJiQkcHByQmppa6TVmZmaws7NTOu7q6qp0zVdffYX9+/fj7NmzGDNmDN5//3388MMPlcayevVqpdVKPT096/7ECKlCSJdu8HN1B1C6JtQXF8OpGUpDDp25g30nogEA+0/ewuWoZ4pzSz8cDGcHa12FRkidUZ8admqd1CxZskRlR93y26NHjzQRq8KXX36J7t27w8/PD5999hkWL16MdevWVVo+NDRUabXSxEQanUI0w5THw6aBQ2FjZo7O7g2wJfgtGp6tAUUlYkTcfIp9J2/h+p0X2PJHhOLcOyM6o2sH3a9BQwjRvlovaLlo0SJMnz69yjLe3t5wc3NDenq60nGpVIrs7Gy4ubmpvM7NzQ1isRi5ublKtTVpaWmVXgMAAQEBWLlyJUQiEczNzSucNzc3V3mcEE3wFNjiwOjxaGbvqJjDhqjXxetPUSKSAgAWfXMI8te1Ye1aemDu+O5VXUqIfmPbjGTkNTW1TmqcnZ3h7Fz9TKeBgYHIzc1FdHQ0/P39AQDnz5+HXC5HQECAymv8/f1hamqK8PBwjBkzBgDw+PFjJCQkIDAwsNLHiomJgb29PSUuRG/4ONJswJoUdqls7pk3CY2lhRmWfzQUJiY8XYVFCGv6vEp3fVDrpKamWrVqhUGDBmHOnDnYtm0bJBIJ5s+fjwkTJsDDwwMAkJycjP79+2P37t3o0qULbG1tMWvWLISEhMDBwQECgQAffvghAgMD0bVr6SiGY8eOIS0tDV27dgWfz8fZs2fxzTff4JNPPtHUUyGE6JGM7AJE36u4wF5RsRjTPt0Nb09H9O/WEm8P8qOmP0KMjMaSGgD4888/MX/+fPTv3x9cLhdjxozBpk2bFOclEgkeP36MoqKyicq+++47RVmRSITg4GBs3bpVcd7U1BRbtmzBwoULwTAMmjVrhg0bNmDOnDmafCqEED1x9nKsonbmvwqKRGjoZo+3+rajhIbUT9T8xAqt0m2Eq5gSUp9N+3Q3nr7MqHCcb26CT2YFYUifNjqIihgyba7S3X7mN+CZsVilW1yCuzuMd5VujdbUEEKIOj1LyFCZ0Hh7OmHlwrfQpKGjDqIihOgLSmoIIfVGWMTDCseG9WuHhTP6gm9uqoOICFEzan5ihZIaQki9IJPLcfpS2RxYlnxTLJ47AAN7tNJhVISoGSU1rFBSQwipF24/SERmTgEAoHkjZ6xc+Ba8PBx0HBUh6kVDutmhpIYQUi+ERZTOTTNqoC8+mtoH5mb09kUIUUbvCoQQvVcikuDmvZdYueAt9O/WUtfhEKI51PzECiU1hBC99yo9D1uWj0dDNztdh0KIRnEYBhwWM62wudYQUFJDCNF73p5Oug6BEFIP0Gp7hJB6IyE1Bzfuv9R1GIRoDqOGzYhRTQ0hRO/J5HLsP30bPx64DL65KfaumQ4HW0tdh0WI2tHoJ3aopoYQotcSUnLw7qr92LjnIkQSGfIKSrDhj391HRYhRA9RTQ0hROdEYimKRRLY2VgAALKFReBxOThx6SG2HbgMkUSmKNva2xUzRwToKlRCNItGP7FCSQ0hRKcYhsE3v57BwEAftGnmjt9PRWHfmVuws7RARnaBopypCQ9zxwRi0uBOMOFRJTMxTNT8xA4lNYQQrcktKIadtYXSsd3HbyLs6iPkFZXg820nUFQsAUcOZJSUJTStvV3x5dxB8G5AC1YSQipHX3cIIVpx+0kSvt55VunYmWuPsHX/ZQDA1TsvUFwoAVcOcF6fN+Fx8cH4Hvh56URKaIhxoNFPrFBNDSFE464/eIlFm/5BW283AKUzBP/49xXsPXVLUea/1eatvd3w5dxgSmaIUaHmJ3YoqSGEaNSlmOdYsuUYxFIZRGIp9p27jV+PXkduTpGiRqY8M1Me5o7phomD/KnvDDE+1FGYFUpqCCEac+7mE3zx00nIZHIAwP3nqbj/LBUcGVQmNADQ068pRvVrTwkNIaTW6F2DEKIRJ68+xP/9eEKR0AAAGIAjrzyhAYDwG0/w9ic7ceTfu5DLjfxrJzFKb5qg6rIZO6qpIYSo3aELd7Fm9zlUWFuvmjdeLocD3xYe6NmxKTr6eILLrSr9IcQAMQwq/uHU8nojRkkNIUSt/jpzCxv+ulDxhLy0lua/LPmm6NquMXp2bIruvk1ga2NRsRAhhNQAJTWEELXZefwGtv59ueIJFc1OAhs+xg/0w9ShnWFmSm9FhAA0+okteichxMCIpVLEZ+QgLi0LXZt5wdFa8ws/MgyDbYevYsex6xXOCSzNISmRwdLKDO1beCAooCX6+DelRIYQVWj0Eyv0rkJIPSWVyZGQlYunqZl4lpaFp2lZiEvNwsusHMjkDPq08sYQ35Yaj4NhGHy/LwJ7zkSjkZs9mns6o4WXM1p4OqO5lzNyhcWwtbaAq6ONxmMhhBg3SmoI0XNyOYPknDxF0hKXloWnaZmIT8+BRCZTeY3AwhzLRvcHh6P5jrY5+cUI6twC80Z1g4W5aYXzLvaUzBBSU5xK+p7V5npjRkkNIXouMi4BP52/juj45BpfY2Vuhi1nI9HA3hYe9gJ42Nuggb0tnG2s1D6iyEFgCQeB5pu4CDEK1PzECiU1hOi57i0aoXuLRniQlIbfL9/CqTtPIJVX/XUsJTcfB2/cr3DchMeFu50NPOwEaGAveJ3wKP9MCCH1FSU1hNQTbRq6Ys2EwQgZ0hN/Xb2DfdfvIq+opFb3kMrkSMzKQ2JWXoVzTjaWuPjFPHWFSwipAxr9xA4lNYToEbmcgVgmg0gqLV0rSSqFWCaDRCqDSCpDK3dnuAis8fGg7pjbrwuO3Y7F75dv43l6ttJ93g/qCmeBFV7lCPEqJx/JOXl4lSNERn5hpXNzedhRLQ0hOkeT77FCSQ0hWlQkluD781dxOe4FxFKZUgIjlsqqbFaa27Mz2jd0U+xbmJliXEB7jO3SDleevMTuy7dw5clLAECGsBAfDAiscA+xVIqU3AKk5AqRnC1Eco6w9OccIXzcndX/hAkhtUI1NexQUkOIFlmamWJJcC8ccXPCt6cvIaeouEbXjfBthYX9u6s8x+Fw0KNlY/Ro2RhxaVn44/JtRDyKh1gqg5kJT6msmYkJGjnZoZGTHdunQggxABEREVi3bh2io6ORkpKCw4cPY+TIkZWWv3z5Mj777DM8evQIRUVFaNSoEebNm4eFCxdqL+gqUFJDiJZxOByM6tAGfVp4Y8PZyzhwq2KH3v/icbk4/fApujbxgp0lv9JyzVwdsXxMEHILiyGSSCskNYQQPafl0U+FhYXw9fXFzJkzMXr06GrLW1lZYf78+Wjfvj2srKxw+fJlzJs3D1ZWVpg7d24dg1YfDsMYXwOcUCiEra0t8vLyIBBQPwKiW9Evk7H8eDiepmdVW5bDAdp6uKKbdyMENvVCR093mJnQdxNCNEkbnxlvHqPr0JUwMa38i0t1pJISXDvxZZ1i5XA41dbUqDJ69GhYWVnh999/r9V1msDV1I2zs7MxefJkCAQC2NnZYdasWSgoKKjymu3bt6NPnz4QCATgcDjIzc1Vy30J0Wf+jRrg0LuT8cmAHuBXs3QAwwD3ktPw06UbmL7rIALW/Ig5vx/GzqvReJyWCSP8jkIIUUEoFCptIpFII49z+/ZtXL16Fb1799bI/WtLY1/xJk+ejJSUFJw9exYSiQQzZszA3LlzsWfPnkqvKSoqwqBBgzBo0CCEhoaq7b6E6DtTHg+ze3TGoDYtsPLkv7j4JF5xblFQD0hkMlx9noCYxBSlzsTFEikuxb3ApbgXAAAna0sEenuhW1MvdPNuBFeBtbafCiGEDTWNfvL09FQ6vGzZMixfvpxFYMoaNmyIjIwMSKVSLF++HLNnz1bbvdnQSPNTbGwsWrdujZs3b6JTp04AgLCwMAwZMgRJSUnw8PCo8voLFy6gb9++yMnJgZ2dndru+wY1PxF9xjAMzsbG4etTF5AmLEBIUHfM7dkFAFAgEiPqZRKuPkvA1WcvEZeRXeW9Dr87Ga3cXbQRNiEGS5vNT4GDv2Ld/BR5aikSExOVYjU3N4e5uXmV19am+Sk+Ph4FBQW4du0alixZgs2bN2PixIl1jltdNFJTExkZCTs7O0XiAQBBQUHgcrm4fv06Ro0apdX7ikQipao3oVBYp8cnRBs4HA4Gtm6Obk0b4Yd/IxGbkqE4Z21uhj4tvNGnhTcAIE1YgMjnpQlO5PMEZBQUKcoK+OZo4eqk9fgJIbonEAg0+qW9SZMmAIB27dohLS0Ny5cvN9ykJjU1FS4uyt8OTUxM4ODggNTUVK3fd/Xq1VixYkWdH5cQXbA2N0PooN5IycuvtIyrwBojO7TGyA6twTAMnqRnIfJZAq4+fwl7SwvwuBrrNkcI0YR6uPaTXC7XWJ+d2qpVUrNkyRKsXbu2yjKxsbGsAtKE0NBQhISEKPaFQmGF9kZC9JW7bc1WueZwOGjp6oSWrk6Y3q2jhqMihGiCtiffKygoQFxcnGI/Pj4eMTExcHBwgJeXF0JDQ5GcnIzdu3cDALZs2QIvLy/4+PgAKJ3nZv369fjoo4/qHrQa1SqpWbRoEaZPn15lGW9vb7i5uSE9PV3puFQqRXZ2Ntzc3Cq5snp1vW9N2hIJIYQQYxMVFYW+ffsq9t9UAEybNg27du1CSkoKEhISFOflcjlCQ0MRHx8PExMTNG3aFGvXrsW8efqxblytkhpnZ2c4O1c/lXpgYCByc3MRHR0Nf39/AMD58+chl8sREBBQt0g1eF9CCCFEL8iZ0o3N9bXQp0+fKqeC2LVrl9L+hx9+iA8//LAukWmFRhrcW7VqhUGDBmHOnDm4ceMGrly5gvnz52PChAmKEUrJycnw8fHBjRs3FNelpqYiJiZGURV27949xMTEIDs7u8b3JYQQQuotRg2bEdNYL8I///wTPj4+6N+/P4YMGYIePXpg+/btivMSiQSPHz9GUVHZaI1t27bBz88Pc+bMAQD06tULfn5+OHr0aI3vSwghhNRXHJT1q6nTpusnoGO0TALNU0MIIaQK2pynpnvQCpiYsJinRlqCK+eWGe3nGy0aQwghhOgLNc0obKwoqSGEEEL0hLaHdBsampmLEEIIIQaBamoIIYQQfVEPZxTWJ5TUEEIIIXqCwzDgsOgXw+ZaQ0DNT4QQQggxCFRTQwghhOgL+euNzfVGjJIaQgghRE9Q8xM71PxECCGkXpHJjLw6glSKkhpCCCF6Le1lhmLRxbzMfEz3WYC/1v6DovxiHUemAbT2EyuU1BBCCNFrx38Ox8lfzgMA/t54AumJWdi1bD++Gr9Rt4FpwpsZhdlsRoz61BBCCNFrL+4n4vAPYfBo6op/fjyjOD7yg2AdRqUZNKMwO5TUEEII0WvP7yVAIpLg82FrIX/dn6aFvzcChvjpODKibyipIYQQorfycwqRmZwNAIqEBgAEjtb4feXf4HA46PpWRzT3a6KrENWLFrRkhfrUEEII0VsvHiSqPB515i7+/OYwXsYmwbudl5aj0hyOnP1mzCipIYQQorcqS2oAoNPA9vhs1wfgmfC0GBHRZ9T8RAghRG/F31ed1LTr4YMv9y6AqZmBfYxR8xMrBvbbQAghxJC8UJHUtPD3xopDi8C3NNdBRBpGq3SzQs1PhBBC9BLDMHjxMEnpWKPWDfH10cWwEljqKCqizyipIYQQDWEYBn/9eRUvX2ToOpR6KSMxC4V5RYp9d29XrD6xBAJHGx1GpVlv1n5isxkzSmoIIUQDxCIp1nx9FL9uv4D/W3IAubmFug6p3okv10nYqYED1pwMhaO7vQ4j0gKaUZgVSmoIIUTNcnML8emiPQg/+wAAkJqSi7/339RxVPXPm/40ts4CrDkVCrfGzjqOiOg76ihMCCFq9PJFBv5vyQGkpuQqjo1+uzOmz+ylu6Dqqef3EmBtZ4nVJ5bAs4WHrsPRDgYAm7lmjLuihmpqCCGkrs6E3VXaj7r5HB++v1uR0HB5HHy0MBjvfzgAPBN6u62ttJcZWHlkMZq2b6TrULSG+tSwQ39lhBBSB7dvvcB360+hpEQCADh6JBqff7YPRYUiAICllTm+XjMew0f66zLMeksul2PW1xPRumtzXYeiXQxY9qnR9RPQLWp+IoSQWiouFmPDupOQSGR4cD8J1yPjcOhgWZ8ZN3c7fL1mLBoZeR8QUbEY96/Hwb9P61pfy+Vy0a6HjwaiIoaMkhpCCKmlXb9GIOVVLgBg1YojyBcWK861adsQK74eAzs7Kx1Fpzs56UKYW5rB0poPANi99hgO/XgO/ccFYN7KsbAxwtek1mhGYVao+YkQQmrh4YNkHDp4Q7FfPqHpP6AN1m2YZJQJDQBcOXEb+zedBgA8vPkMh7eFAwDC919HVPgDXYZWf8jVsBkxqqkhhJAaEoul+N+3J1R+GW7v64lFi4fCzNDWIqqFyLA7uBf5FP3e7oINH+0G8/qFChjYDn1Gd9ZxdMQYUE0NIYTU0J7fr+Dli0yV5+7eScQ747fgt50RyCs3C66xKBQW4+6VJ5CIpFj01jokP08HAFjbWeKj9ZPB4XB0HGH9QKOf2DHerxSEEFILz5+l468/I6ss09DLEc2aucL6dZ8SYxJ1/gGkEhkAoCCvrEnuva/HwcHVVldh1T/Up4YVSmoIIaQaMqkc69ceh0xWscOCqSkP/fq3wai3O6NZc1cdRKcfroXdUXl80yd7sG/TaTRt64kP1k6AlY2FliMjxoSSGkIIqcbBAzfw5HGq0jEHBysMG9ERbw33g72DtY4i0w9SiQw3K+kILCoWQyaRYVLIYEpoaoJqaljRWJ+a7OxsTJ48GQKBAHZ2dpg1axYKCgqqvGb79u3o06cPBAIBOBwOcnNzK5Rp3LgxOByO0rZmzRoNPQtCiLFLSszCbzsiFPstWrrhs8+H4Y99H2DK9J5Gn9AAwL2rT1BYbhRYeX69W+G7U4vRsJmblqOqp2hBS1Y0VlMzefJkpKSk4OzZs5BIJJgxYwbmzp2LPXv2VHpNUVERBg0ahEGDBiE0NLTScl999RXmzJmj2LexMdxl6AkhuiOXM/jftychlcrQq7cPRo/tjDZtG1Kn1/+I/M9yEW+MnNsPs5eNBs+Ep+WIiLHSSFITGxuLsLAw3Lx5E506dQIA/PDDDxgyZAjWr18PDw/VC5MtWLAAAHDhwoUq729jYwM3N8r6CSGadfHfh2jVugGW/N9wuLpRZ1dVGIbBtdPK/WlMTHmY/+0kBE/qpqOo6jE5ADY5s5HPU6OR5qfIyEjY2dkpEhoACAoKApfLxfXr11nff82aNXB0dISfnx/WrVsHqVRaZXmRSAShUKi0EUKqViyW4NCVewiLeqzrUHSmV59WmPteP0poqvD8QRIyknMU+7aO1lj99wJKaOqIhnSzo5GamtTUVLi4uCg/kIkJHBwckJqaWslVNfPRRx+hY8eOcHBwwNWrVxEaGoqUlBRs2LCh0mtWr16NFStWsHpcQoxFUmYu9kfcxT+R9yEsEuHI0uk1vragWIQSsRROtoYxoy6PR1N5Vaf8qCfvNg2x9Ld34erpqMOI6jnqKMxKrZKaJUuWYO3atVWWiY2NZRVQdUJCQhQ/t2/fHmZmZpg3bx5Wr14Nc3NzldeEhoYqXScUCuHp6anROAmpT+RyBpGxL7EvIgaXH8Qr3hfbN3FHI1d7yOUMcguLkZ5bgLTcfKTlFCAjrwBpOQVIzytAek7p8SKRBD3aNMEP74/U6fMh2vOmP033oR2w6IdpsLAyvjl6iP6oVVKzaNEiTJ8+vcoy3t7ecHNzQ3p6utJxqVSK7OxstfeFCQgIgFQqxYsXL9CyZUuVZczNzStNeAgxZvlFJTh67SH2RdxBYkZuhfPJmXl4a+mvSM8rhEQqq9E903Lz1Rwl0VfpSdl4di8Rkz8ZikmLhoDLpZot1uQMwGFR2yKnmpoac3Z2hrOzc7XlAgMDkZubi+joaPj7+wMAzp8/D7lcjoCAgLpFWomYmBhwudwKzV2EkMo9Tc7EvogYnLgRixJx5X3SsvJrN92/tYU5rMzN2IZH6omYS4/w+c+z0XO4v65DMRzU/MSKRvrUtGrVCoMGDcKcOXOwbds2SCQSzJ8/HxMmTFCMfEpOTkb//v2xe/dudOnSBUBpX5zU1FTExcUBAO7duwcbGxt4eXnBwcEBkZGRuH79Ovr27QsbGxtERkZi4cKFeOedd2Bvb6+Jp0KIwZDIZLhw5xn2RdxB9NOkWl/vYGMJFztruNhZw/X1v6U/28DV3houttaw5FNCY0x6jegEviX9P6/PIiIisG7dOkRHRyMlJQWHDx/GyJEjKy1/6NAh/Pjjj4iJiYFIJEKbNm2wfPlyBAcHay/oKmhsnpo///wT8+fPR//+/cHlcjFmzBhs2rRJcV4ikeDx48coKir7Jrht2zalDr29evUCAOzcuRPTp0+Hubk59u7di+XLl0MkEqFJkyZYuHChUn8ZQohqmXmFKBJJ0LJhaW3rk6QM5BeLqrzGhMfFrwvHwaehM8xMaQJyoowSGk1gO4Fe7a4tLCyEr68vZs6cidGjR1dbPiIiAgMGDMA333wDOzs77Ny5E8OGDcP169fh5+dX16DVhsMwxldXJRQKYWtri7y8PAgEAl2HQ4hOMAyDlOx8PEnOwOOkDDxJzsCTpAwkZeYplQv2b4k1M4foKEpCdE8bnxlvHiOoyYcw4da9D6hULsK5+B/qFCuHw6m2pkaVNm3aYPz48Vi6dGmtrtME+upFiJHicDjwcBTAw1GAPu2bKo4XFIvwNDlTkeg8TsrA9UcJCPDx0mG0hJDa+O98bJoaMCOXy5Gfnw8HBwe137suKKkhhCixtjCHX7MG8GvWQNehEGJ85Axq24RU8XpUmLZk2bJlWL58ed3vW4n169ejoKAA48aNU/u964KSGkIIIURfMPLSjc31ABITE5WanzRRS7Nnzx6sWLEC//zzj96MQKakhhBCCDEwAoFAo31G9+7di9mzZ+PAgQMICgrS2OPUFiU1hBBCiL6oB/PU/PXXX5g5cyb27t2LoUOHavzxaoOSGkIIIURfqKlPTU0VFBQo5oYDgPj4eMTExMDBwQFeXl4IDQ1FcnIydu/eDaC0yWnatGn4/vvvERAQoFjP0cLCAra2ul/4lea0JoQQQvTFm5oaNlstREVFwc/PTzHHTEhICPz8/BTDs1NSUpCQkKAov337dkilUnzwwQdwd3dXbB9//LH6XgMWqKaGEEIIMVJ9+vRBVdPV7dq1S2n/woULmg2IJUpqCCGEEH3BgGWfGrVFUi9RUkMIIYToi3rQUVifUZ8aQgghhBgEqqkhhBBC9IVcDoDF5HtyFtcaAEpqCCGEEH1BzU+sUPMTIYQQQgwC1dQQQggh+oJqalihpIYQQgjRF1qeUdjQUPMTIYQQjcvPKdR1CMQIUFJDCCFE475f+DsykrMBAPEPk1BSKNJxRPqJYeSsN2NGSQ0hhBCNSk/KwtXjt/DP9vNIT8rCkpEbsHDQGqS+zNB1aPqHYUqbkOq6GXmfGkpqCCGEaNSp3y5BLmdwclcEVryzFXmZ+Yh/kIQV72yF3MjnValAywtaGhpKagghhGiMRCxF2B+XAQBF+cV4drd0xWdTMxN8/N0UcLn0MUTUh0Y/EUII0ZjIE7eRk5ZX4XjrgKYQiyQoKRSBb2Wug8j0lFwOcFjUXhl5nxpKagghhGjM8Z0XVR6/c+kx7lx6DCd3O6w6sACNWzfQcmR6imE5pNvIm58oqSGEEKIRCY9TcPfy40rP9xjujw//Nxm2jjZajIoYMkpqCCGEaMSJSmpprO0s8cG6Segzugs4HI6Wo9JvjFwOhkXzk7EP6aakhhBCiNqVFIpwbu/VCsc7B7XFgu+nwdHdTvtB1QfU/MQKJTWEEELU7sKhGygUFiv2LazNMXfVeAya0oNqZ4jGUFJDCCFErRiGwfEdFxT77bq3wKLN0+HWyFl3QdUXcgbgUE1NXVFSQwghRK2e3HqBuDsJMOObYsaXozFiXj+aj6amGAYAmyHdlNQQQgghanN8xwW07NgEn2ydAc8W7roOhxgRSmoIIYSoTXFBCRq3boAF308Fz4Sn63DqHUbOgGHR/MRQTQ0hhBCiHhbWfIz5YKCuw6i/GDnYNT/RkG5CCCGE6AGqqWGHem4RQgghxCAYZU3Nm0xWKBTqOBJCCCH67s1nhTZqQaSMiFUTkhQSNUZT/xhlUpOfnw8A8PT01HEkhBBC6ov8/HzY2tpq5N5mZmZwc3PD5dSTrO/l5uYGMzMzNURV/3AYI2yAk8vlePXqFWxsbCqd2VIoFMLT0xOJiYkQCARajrB+odeq5ui1qjl6rWqOXqvaqe3rxTAM8vPz4eHhodH5dkpKSiAWi1nfx8zMDHw+Xw0R1T9GWVPD5XLRsGHDGpUVCAT0JlFD9FrVHL1WNUevVc3Ra1U7tXm9NFVDUx6fzzfaZERdqKMwIYQQQgwCJTWEEEIIMQiU1FTC3Nwcy5Ytg7m5ua5D0Xv0WtUcvVY1R69VzdFrVTv0ehkuo+woTAghhBDDQzU1hBBCCDEIlNQQQgghxCBQUkMIIYQQg0BJDSGEEEIMAiU1hBBCCDEIRp3UbNmyBY0bNwafz0dAQABu3LhRZfmNGzeiZcuWsLCwgKenJxYuXIiSkhItRatbtXmtJBIJvvrqKzRt2hR8Ph++vr4ICwvTYrS6ExERgWHDhsHDwwMcDgdHjhyp9poLFy6gY8eOMDc3R7NmzbBr1y6Nx6kPavtapaSkYNKkSWjRogW4XC4WLFiglTj1QW1fq0OHDmHAgAFwdnaGQCBAYGAgTp8+rZ1gday2r9Xly5fRvXt3ODo6wsLCAj4+Pvjuu++0EyxRO6NNavbt24eQkBAsW7YMt27dgq+vL4KDg5Genq6y/J49e7BkyRIsW7YMsbGx+PXXX7Fv3z58/vnnWo5c+2r7Wn3xxRf46aef8MMPP+Dhw4d49913MWrUKNy+fVvLkWtfYWEhfH19sWXLlhqVj4+Px9ChQ9G3b1/ExMRgwYIFmD17tlF8ANX2tRKJRHB2dsYXX3wBX19fDUenX2r7WkVERGDAgAE4efIkoqOj0bdvXwwbNoz+BlWwsrLC/PnzERERgdjYWHzxxRf44osvsH37dg1HSjSCMVJdunRhPvjgA8W+TCZjPDw8mNWrV6ss/8EHHzD9+vVTOhYSEsJ0795do3Hqg9q+Vu7u7szmzZuVjo0ePZqZPHmyRuPUNwCYw4cPV1lm8eLFTJs2bZSOjR8/ngkODtZgZPqnJq9Veb1792Y+/vhjjcWjz2r7Wr3RunVrZsWKFeoPSI/V9bUaNWoU884776g/IKJxRllTIxaLER0djaCgIMUxLpeLoKAgREZGqrymW7duiI6OVjS7PH/+HCdPnsSQIUO0ErOu1OW1EolEFRZls7CwwOXLlzUaa30UGRmp9NoCQHBwcKWvLSF1IZfLkZ+fDwcHB12Hovdu376Nq1evonfv3roOhdSBUa7SnZmZCZlMBldXV6Xjrq6uePTokcprJk2ahMzMTPTo0QMMw0AqleLdd981+OanurxWwcHB2LBhA3r16oWmTZsiPDwchw4dgkwm00bI9UpqaqrK11YoFKK4uBgWFhY6iowYkvXr16OgoADjxo3TdSh6q2HDhsjIyIBUKsXy5csxe/ZsXYdE6sAoa2rq4sKFC/jmm2+wdetW3Lp1C4cOHcKJEyewcuVKXYemd77//ns0b94cPj4+MDMzw/z58zFjxgxwufTrRoi27dmzBytWrMD+/fvh4uKi63D01qVLlxAVFYVt27Zh48aN+Ouvv3QdEqkDo6ypcXJyAo/HQ1pamtLxtLQ0uLm5qbzmyy+/xJQpUxTZe7t27VBYWIi5c+fi//7v/wz2A7sur5WzszOOHDmCkpISZGVlwcPDA0uWLIG3t7c2Qq5X3NzcVL62AoGAamkIa3v37sXs2bNx4MCBCs2cRFmTJk0AlL63p6WlYfny5Zg4caKOoyK1ZZifxNUwMzODv78/wsPDFcfkcjnCw8MRGBio8pqioqIKiQuPxwMAMAa8JmhdXqs3+Hw+GjRoAKlUir///hsjRozQdLj1TmBgoNJrCwBnz56t9rUlpDp//fUXZsyYgb/++gtDhw7VdTj1ilwuh0gk0nUYpA6MsqYGAEJCQjBt2jR06tQJXbp0wcaNG1FYWIgZM2YAAKZOnYoGDRpg9erVAIBhw4Zhw4YN8PPzQ0BAAOLi4vDll19i2LBhiuTGUNX2tbp+/TqSk5PRoUMHJCcnY/ny5ZDL5Vi8eLEun4ZWFBQUIC4uTrEfHx+PmJgYODg4wMvLC6GhoUhOTsbu3bsBAO+++y42b96MxYsXY+bMmTh//jz279+PEydO6OopaE1tXysAiImJUVybkZGBmJgYmJmZoXXr1toOX6tq+1rt2bMH06ZNw/fff4+AgACkpqYCKO2wb2trq5PnoC21fa22bNkCLy8v+Pj4ACgdDr9+/Xp89NFHOomfsKTr4Ve69MMPPzBeXl6MmZkZ06VLF+batWuKc71792amTZum2JdIJMzy5cuZpk2bMnw+n/H09GTef/99JicnR/uB60BtXqsLFy4wrVq1YszNzRlHR0dmypQpTHJysg6i1r5///2XAVBhe/P6TJs2jendu3eFazp06MCYmZkx3t7ezM6dO7Uety7U5bVSVb5Ro0Zaj13bavta9e7du8ryhqy2r9WmTZuYNm3aMJaWloxAIGD8/PyYrVu3MjKZTDdPgLDCYRgDbjshhBBCiNEwyj41hBBCCDE8lNQQQgghxCBQUkMIIYQQg0BJDSGEEEIMAiU1hBBCCDEIlNQQQgghxCBQUkMIIYQQg0BJDSGEEEIMAiU1hBBCCDEIlNQQQgghxCBQUkMIIYQQg/D/dJtevnK8QKgAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "r, z = np.random.multivariate_normal(\n", " mean = [1.1, 0], cov = [[0.01, 0], [0, 0.01]],\n", " size = 50\n", ").T\n", "\n", "x = r\n", "y = 0 * r\n", "\n", "Bx, By, Bz = field.interpolateXyz([x, y, z])\n", "\n", "plt.quiver(r, z, Bx, Bz, By)\n", "plt.colorbar(label = \"$B_T$\")" ] }, { "cell_type": "code", "execution_count": null, "id": "e45c70a7-5020-4877-89be-b19b1152aa3c", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.4" } }, "nbformat": 4, "nbformat_minor": 5 }